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SECTION I 
INTRODUCTION 


The F0RTRAN IV program described herein analyzes the nozzle for a hypersonic 
scramjet by a two dimensional second order characteristic procedure described 
In ATL TR 186, “A Source Flow Characteristic Technique for the Analysis of 
Scramjet Exhaust Flow Fields.'* 

The program starts from the Initial profile and marches along down-running char- 
acteristics until the final X station Is reached. This process can be interrupted 
by "Change of Origin Profiles" in source flow cases. This is done by interpolat- 
ing the characteristic data at the required axial station and setting up "F 
a-rays" i.e., XF, YF, ... . Then the program continues on down-running char- 
acteristics using the "F arrays" as an initial profile. These "F arrays" are 
computed for all flows at an X station that coincides with the cowl tip if 
there is a cowl in the flow field. If the flow is overexpanded at the cowl the 
program will invert the problem as described in ATL TR 186. When the shock 
reflects off the lower wall "F arrays" are again calculated and the computation 
continues with the problem inverted again until the final axial (X) station is 
attained. 
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SECTION II 

DESCRIPTION OF INPUT 


A, Input Format 


Card 

Number Columns Format Description 


1 


1-5 


15 


Jl, type of flow (0**two dimensional, 1-axisym- 
metric 


6-10 

n-15 

16-20 

21-25 

26-30 


15 J2, coordinate exponent for line source system 

(O-two dimensional or axisymmetric, 1-line 
source) 

15 NPTS, number of data points on Initial profile 

15 lEQ., chemistry indicator (0-frozen hydrogen- 

air chemistry, 1 -equi 1 tbrium chemistry) NOTE: 
if IC0WL “ 1 then I EQ = 1 . 

15 IC0WL, external data Indicator (1 -overexpans ion 

or underexpansion Interaction calculations Is 
required, O-Internal flow only)* If lEQ »= 0 
then IC0WL = 0. 

15 I0VER, overexpansion indicator (O-flow de- 

finitely underexpanded, 1-flov; overexpanded 
or marginal) 


2 


1-10 

ElO.O 

•'XBP , ratio of axial coordinate of cowl at Ini- 
tial station to throat height 

11-20 

ElO.O 

XB0D, ratio of axial coordinate of lower wall 
at initial station to throat height 

21-30 

EiO.O 

XC0WL, ratio of axial coordinate of cowl trail- 
ing edge to throat height 

31-^0 

E10*0 

RTH, throat height (ft.), scaling parameter L* 

^1-50 

EIO.O 

TEST, maximum allowable axial step size, used 
for computing upper boundary - typical value 

= .1 


*(A11 length variables are nondi mensicnal ized by RTH,) 


2 - 
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Card 


Number 

Columns 

Format 

Description 

2 

...51-60 

ElO.O 

XFINAL, ratio of final axial coordinate of run to 
throat height 


61-70 

ElO.O 

XTJl, ratio of axial coordinate of beginning 
of Cartesian region to throat height 

3 

1-10 

ElO.O 

XSHFT, ratio of axial coordinate of moment 




axis to throat height 


11-20 

ElO.O 

YSHFT, ratio of radial coordinate of moment 
axis to throat height 


.21-30 

ElO.O 

XTHX, initial thrust (Ibs/RTH^) 


31-40 

ElO.O 

YLFT, initial lift (Ibs/RTH^) 


1-5 

15 

NXXJl, number of locations in source flow 
where a new initial profile Is desired, maxi- 
mum is 4, minimum is 1 (I.e. changes of origin, 
cowl station,,..) 

5 

One of the 

following 

cards Is necessary for each NXXJl. If NXXJl equals 


1 and 2-D 

flow, card 

5a may be blank. 

5a 

1-10 

ElO.O 

XXJl, ratio of axial coordinate of new Initial 
profile to throat height 


11-20 

ElO.O 

AXX, coefficients of polynomial describing a 
segment of lateral extent of the nozzle 


21-30 

ElO.O 

BXX, for the equation Z = AXX(X-X,) + BXX 

^ * 


31-40 

ElO.O 

XOR, ratio of axial coordinate of line source 
origin 

5b 

1-5^ 

15 

1 FENCE, fence indicator (0-no fence;) 
1-supersonic fence exists) 


11-20 

ElO.O 

AFENCE, coefficients of fence for the equation 


21-30 

ElO.O 

BFENCE, Y=AFENCE (X-XBP) + BFENCE 


31-40 

ElO.O 

XFENCE, ratio of axial coordinate of fence on 
lower surface to throat radius 

6 

1-5 

15 

NUWSEG, number of polynomial segments describ- 
ing the cowl (maximum is 3) 


6-10 

15 

NLWSEG, number of ploynomlal segments describ- 
ing the lower wall (maximum Is 3)* 




- 3 - 
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Card 


'Number 

Col umns 

Format 

Description 

7 

One of 

the following 

cards Is necessary for each cowl segment. 

7A 

1-10 

ElO.O 

XXU(L), ratio of axial coordinate of begln- 
ing of segment of cowl to throat height 


11-20 

E10*0 

t h 

A, coefficients of the “L “ segment 


21-30 

E10*0 

B, of polynomial describing cowl for the 




- . equation 


31-40 

ElO.O 

C, Y«AX^4-BX+C; Y=Y/RTH S X^X/RTH 

8 

One of 

the following 

cards. is necessary for each lower wall segment 

8a 

1-10 

ElO.O 

XXL(L), ratio of axial coordinate of beglnlng 

of segment-- of lower wall to throat 

height 


11-20 

ElO.O 

t h 

A, coefficients of the "L'' " segment 


21-30 

ElO.O 

B, of polynomial describing lower wall for 




the equation 


31-40 

ElO.O 

C. Y-AX^+BX+C; Y=Y/RTH S X=X/RTH 

. 9 

1-10 

ElO.O 

EMINF, free stream or reference Mach number 


11-20 

ElO.O 

TIN, free stream or reference temperature ( K) 

/ 

21-30 

ElO.O 

WINF, free stream or reference molecular 
Weight (lbs/1 b-mole) 


31-40 

ElO.O 

PlNF, free stream or reference pressure 
(Ibs/in^) 

10 

Initial Profile - One of the following cards are required for each 
data point as described below* The program reads data points be- 
ginning at the lower wall (point #1) and proceeds to the cowl 
(point #NPTS). 

10a 

1-10 

ElO.O 

X(l), ratio of axial coordinate of data 
point to throat height 


11-20 

ElO.O 

Y(l), ratio of radial coordinate of data 
point to throat height 


21-30 

ElO.O 

P(l), ratio of pressure at data point to 
free stream pressure 


31-40 

ElO.O 

0,(1), ratio of velocity at data point to 


free stream velocity 

-ii- 



Col urns 

Format 

Descri ption 


ElO.O 

* 

T(I), ratio of temperature at data point to 
^ free stream temperature 

51-60 

ElO.O 

TH(I), flow Inclination at data point (In 
rad 1 ans) 

61-70 

ElO.O 

R(l), fuel to air equivalence ratio (only 
necessary If lE(i=l) 

Initial Mass Fractions - Same as card 10 above except for mass 
fractions at each data point. These cards are necessary only 
if IEQ-0. 

1-10 

ElO.O 

mass fraction of H 

11-20 

ElO.O 

mass fraction of 0 

0 

1 

CM 

ElO.O 

mass fraction of H^O 

31-^0 

ElO.O 

mass fraction of H 2 

i»l-50 

ElO.O 

mass fraction of 0^ 

51-60 

ElO.O 

mass fraction of OH 

61-70 

ElO.O 

mass fraction of N 2 

This card Is necessary only If 1C0WL=I. 

1-5 

15 

MM, number of points in Prandtl -Meyer fan 
including data downstream of shock(MM~9) 

11-20 . 

ElO.O 

PM, ratio of pressure external to the cowl 
to free stream or reference pressure 

21-30 

ElO.O 

QM, ratio of velocity external to the cowl 
to free stream or reference velocity 

' 31-AO 

ElO.O 

TM, ratio of temperature external to the 
cowl to free stream or reference temperature 

iil-50 

ElO.O 

THM, flow inclination external to the cowl 
(In radians) 

51-60 

ElO.O 

RM, fuel to air ratio external to the cowl 
(only needed If IEQ»1) 




FIGURE #2« Definition of Change of Origin Variables 

- 6 - 
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SECTION 111 

DESCRIPTION OF OUTPUT 

A* Output Format - The heading page contains program constants, line source 
coordinates and origin changes, upper and lower wall coordinates, free stream 
data at the cowl and Infinity conditions. The program then prints the "Initial 
profile," the data at selected points along down-running characteristics In the 
flow field and the running lift, thrust and pitching moment. The process con- 
tinues until the input value for the final axial location Is reached. 

The above flow can be interrupted by "change of origin profiles" and a profile 
at the trailing edge of the cowl. If the flow Is ovei expanded at the cowl the 
output will switch to up-running characteristics from the lower wall to the con- 
tact surface and additional output covering the shock angle and external Mach 
number will be printed. When the shock reflects off the lower wall another pro- 
file will be printed and the run will proceed to the final axial location through 
down-running characteristics. 

B. IdentI fi elation of Output VcTriables 

X - axial distance/throat height 
Y - radial distance/throat height 
Q - veloci ty/f ree stream velocity . 

T - temperature/f ree stream temperature 


p 

- 

pressure/f ree 

stream 

TH 

.. 

flow 

angle (radians) 

EM 

- 

Mach 

number 



ALP(l) 

= 

mass 

fraction 

of 

H 

ALP (2) 

s: 

mass 

fraction 

of 

0 

ALP (3) 


mass 

fraction 

of 


ALP(1() 

rs 

mass 

fraction 

of 

«2 

ALP (5) 

= 

mass 

fraction 

of 

°2 

ALP (6) 

= 

mass 

fraction 

of 

OH 

ALP (7) 

= 

mass 

f ract.i on 

of 

«2 


- 7 “ 
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Frozen flow extra output Is: 

CPX - specific heat/free stream specific heat 
W - molecular weight of mlxture/free stream molecular weight 

Equilibrium flow extra output Is: 

GAM - equilibrium Isentropic exponent 
PHI - fuel to air ratio 
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SECTION IV 

SUBROUTINES AND FUNCTIONS 


A. Subroutines 

Name Description 


K INDATA 

2. COEFF 

3. ERROR 

4. BODL 

5. BODU 
COWL 

7. COWLO 

8. SHOCK 

9. SHOCPT 

10. PANDW 

1 1 . DRTEST 

12. PM 

13* all 

14. THERMO 

15. XMASSS 

16. FM 

17. PMl 


reads and prints initial data and computes infinity conditions 

sets theromodynamic coefficients as functions of temperature 

prints program statement number nearest selected errors and 
terminates computer run 

locates axial and radial coordinates along lower vehicle surface 

locates axial and radial coordinates along upper cowl surface 

calculates shock jump relations and Prandtl~Meyer fan at cowl 
trailing edge for under-expanded flow 

computes contact points, shock points and shock angle at cowl 
trailing edge for over-expanded flow 

calculates shock jump conditions 

computes shock angle by matching pressure from shock jump re- 
lations and pressure from characteristic relation on downstream 
side of shock 

makes pressure and frozen chemistry molecular weight or equili- 
brium chemistry fuel to aiV ratio 

tests for dropping data points on free stream side of shock 
computes i sen tropic air expansion 

calculates density, ratio of specific heats, Mach number, Mach 
angle, frozen flow specific heat, molecular weight and gas 
constant 

computes frozen flow specific heat, derivative of specific heat 
and enthalpy of each species from polynomial fits in temperature 

computes mass flow correction factor 

computes mass function 

computes flow properties for given Prandtl -Meyer turning 


- 9 - 
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B. Functions 


Name 
1 » FT 

2. FH 

3. FGAM 

RHEQ 

5. XMl 

6. XM2 
7* XM3 
8. FI 

3. F2 


/ 


Description 

calculates temperature at any data point in flow field 

computes static enthalpy at any point In flow field 

makes Isentropic exponent and associated thermodynamic data at 
any point In flow field 

computes equilibrium density 

calculates tan (6+p) along up-running characteristic 

calculates tan (6-y) along down-running characteristic 

calculates tan (0) along streamline 

or Bj coefficient along up-running and down-running 
characteristics- respectively 

*A^ or coefficient along up- running and down- running 
characteristics respectively. 


*Note: p - e relationship along up-running characteristic: 

A,(Pc-P;^) + ec'V^a^^c'V “ ° 

p - 8 relationship along down-running characteristic: 


- 10 - 
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MACHINE CONTROL CONSIDERATIONS 

Machine - program designed for CDC 6600. 

Estimates for run. 

a. Field length: 

(1) compile - 65,000 octal locations 

(2) load - 120,000 octal locations 

(3) run - 100,000 octal locations 

b. CP time: variable depending on number of points in Initial 
profile and type of flow 

c. 10: less than 100 octal seconds 

d. Tapes or disks used: 

(1) Tape 5 ” card input 

(2) Tape 6 - printed output 

(3) no other tapes or disk files used 

e. Printed output: as In CP time It is variable depending on 
case submitted. 






TR 179 • 

. 


1 

- ' - 

^ ■ 

♦ 

./... SECTION . IV 



. , 

— 

SAMPLE . 

INPUT FOR 

SOURCE FLOW CASE 

IB FROM ATL 

TR 186 -- 



0 

1 4 1 

1 1 

0 




0 • 

0 , 

3 . 

1 . 

.1 

21 . 

22 . 

0 . 

.5 

0 . 

vO - 

0 . 


• 

2 







0 . 

,0 875 

1 . 

-7 . 




3 , 

.0875 

1 . 

-7 . 




1 

59 

3 . 064 

16 . 




2 

3 






a . 

.1314 

0 . 

1 * 




0 . 4 

0 . 

.1051 

.979 




0 . 

-.55 65 

0 . 

0 , 




0 . 4 

0 , 

4452 

.0 89 05 




8 . 

.01019 

60 82 

.7 41 




10 . 

232.3 

28.9 6 

.1603 




0 . 

0 . 

36 .6 5 

.929 

10 .1 

0 . 

1 - 

0* 

.0 25 

36 .6 5 " 

.9 29 

10 .1 

0. 

1 . 

0. 

.0 5 

36 .6 5 

.929 

10 .1 

0 . 

i * 

0 , 

.0 75 

36 .65 

.9 29 

10 .1 

0 . 

1 . 

0 . 

.1 

36 .65 

.929 

10 .1 

0. 

1 . 

0 * 

.125 

36 ,65 

,929 * 

10 .1 

0 . 

1 . 

0, 

.15 

36 .6 5 

.9 29 

10 .1 

0 , 

1 . 

0 . 

.175 

36 .6 5 

.929 

10 .1 

0 . 

1 . 

0. 

.2 

36 .6 5 

.9 29 

10 .1 

0. 

1 . 

0 , . 

.2 25 

36 .6 5 

.9 29 

10 .1 

0 . 

1 . 

0 . 

,25 

36.65 

.9 29 

10 .1 

0 . 

1 . 

0 . 

,2 75 

36 .6 5 

.929 

10 ,1 

0 . 

1 . 

0 , 

.3 

36 .65 

.929 

10 .1 

0 , 

1 , 

0 . 

.325 

36.65 

.929 

10.1 

0. 

1 . 

0 . 

.35 

36.65 

.929 

10.1 

0 . 

1 . 

0 . 

.375 

36,65 

.929 - 

10.1 

0 . 

1 . 

0 . 

.4 

36 .65 

.929 

10 .1 

0 . 

1 . 

0. 

.4 25 

36 .65 

,929 

10 .1 

0 . 

1 . 

0. 

.4 5 

36 .6 5 

.9 29 

10 .1 

0 . 

i . 

0 . 

.4 75 

- 36.6 5 

,9 29 

10 ,1 

0 . 

1 . 

0 . 

.5 

36 .6 5 

.9 29 

10 .1 

0 . 

1 . 

0 . 

.525 

36 .65 

.929 

10 .1 

0 . 

1 . 

0. 

.55 

36 .65 

.929 

10 .1 

0 . 

1 * 

0. 

.575 

36 .65 

.929 

10 .1 

0. 

1 . 

0 . 

. .6 

36 .65 

.9 29 

10 .1 

0 , 

1 . 

.0* 

.625 

36,6 5 

.929 

10 .1 

0 . 

1 . 

6« 

.6 5 

36 , 6 5 

,929 

10 .1 

0 . 

~ i; ■' 

0 . 

.675 

36 .6 5 

.929 

10 »1 

0 . 

1 . 

0 . 

.7 

36 .65 

.9 29 

10 .1 


1 . 

A * 

.725 ■ 

36.65 

,929 

10 .1 

0 . 

1 . 

0 . 

.75 

36 .6 5 

.929 

10 .1 

0 . 

1 , 

0 . 

.775 

36.65 

.929 

10.1 

0 . 

1 . 

0 * 

.8 

36 ,65 

• .929 " 

10 .1 

0 . 

'i: 

0 , 

.8 25 

* 36.65 

.929 

10.1 

0 . 

1 . 

0 . 

.85 

36 . 6 5 

.9 29 

"10 ,1 

0 . ■ ■ 

'■■" 1 , 

_pj 

^,8 75 

36 .65 

.929 

10 .1 

0 . 

1 . 

0 . 

’ ,9 

36 .6 5 

.929 

10 .1 

0 . 

1 . 

0 , 

.925 

36 .65 

.929 

10.1 

0 . 

1 . 

0 « 

.9 5 

36 .6 5 

.9 29 

10 :i 

■' 0 . 

1 ^ 

0 . 

,975 

36.65 

.9 29 

10,1 

0 . 

1 . 

0 . • 

1 . 

36 .6 5 “ 

.929 

io .1 

0 . 

1. 

9 . 

1 . 

1 . 

A » 

0 , 

0 . . 
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PROGRAM LISTING 

* • \ *- 

The following is a listing of the Fortran IV program for Source Flow 
Character 1st ic 




PROGRAM NOZ BOD ( I NPUT , OUTPUT, PU_NCH^TAPE5= IN_PuTj2A^6=OUTPUT_j-_ _ 

tAPE7 = PUNCH'r" " "■ ■■ • 

CO_MMON/COWL/_IC_0_WLf MM ,XM< 9) »YM< P)_r_PMM (9 )_, WM (9 ) , RHM(9) ,THM( 9), QM_( 9 ) , 
■”YrM fPT, TM <9TT GM (9T; WQM (9y7Er4M( 9T , ALP~Mr779 

COMMON/1 PP /I PP 

COMMON /SHAPE/A AT'{372 ) , AA'2( 3.2) , A A3T3 V2 7TX U< 3) ,X X'l (3 ) , X I NT U, X I NT L 

CO MM ON /L TM /X SH FT ,Y $H FT ,X JH X. YL FT ,X MO M ___ __ 

COMMON/ X'X J /N X X J 1 , X X J 1 ( 6 ) , A X X ( 6 ) f B X X ( 6 > , X ORI 6’ )‘“ 

qOMMON/XFl NA_L/ XP inal_ '■ 

COMMON /rEQ'7l EQTP INVR HFINF/uT NF ,P"INF 

COMMQN/V AR /RpO (2 00 )» ' _ _ 

£EM (2 0a V7XMUr20 0) ,CPX"(200')', W( 20 0) .R (2 00 ) , GAM( 20 0) ,XMASS (2 00"). 

2 _ XN (2 00 )f YN (2 00 ). QN (2 00 ). TN (2Q0 >* PN (2 00 ). THN( 20 0) »RHON 

3 (2W")r'EM N(’20 0) ,XMUN( 20 0) ,CPXN( 20 0) , W N( 20 0) ,R N( 20 0)‘,G”AM N( 20 0) , 

4 AS SN_( 2 00 AL_^_a Oj- 20 0 ) ,S 1 (_102j» Hm 0 )_j_ ER R ( 20 ). T EMP ( 2 0) 

57A'LP'DUMrioy 

COMMON /SP/ _ 

CO MM ON 7F VA R/' ' ' 

IRHOP (2 00 ). CPXF (2 00 ). EMF( 20 OJ ,XMUF( 20 0) ,WF( 20 0) >RP( 2 0 OJ ,GAMF( 20 0) , 

■ 2 XM AS SF (2 00 ) , ""H F ( 20 0) ,S F ( 20 O') , A LP DH-O )7' 

3 tH ET A( 20 ) ^ _ _ , ^ 

COMMON 7X7~XT20'0r;Y (2 00T‘. P'( 2'C orTO (2 OOl. T( 20 0) ,tH( 20 "0) ,ALP (1 d. 20 0) 

common / 1/ lOPuT 

CO MM ON /A / T 1 N , CP I N ,'R"0" 

COMMON/8/ WTMOLE 

COMMON/D/ 'GAM I NF , EM I NF 7R I NF ,'w I NF 

C0MM0N/F/A9,B9. I80D,XWF,NBOD,YEND __ 

' C0MM0¥/6 tX/X J , X Ji yNP'tS n 0 .'I'R'EF E , Tc'hE'M yTPlJNCH'V I DE SON .TRTnxX , X'bP . 
lYBP.THeP,RAD,XBOD, Y80D.THBOD,RAD8,XEND,THEiliyt^TH,YEXI T ,THST.TEST 
' 1 1 R FL , Y 0 , R A 08 2 .'R R A D' (2 0 ) 7 N S T A R . Y N 0 Z E I N P E N , H 1 6 , fU 7 
_ COMMON/XF/XF (20 0 ) , YF ( 20j0 ) , P F ( 20 0 )_. OF (200 ) , T_F (JOO ) , T_H F ( 2 00 ) , ALPF (_ 

‘““iio‘.2oo ) ~ ' * 

_cO mm oN /X Cp WL /X cp>i 

' CO mm ON /C'T/ ‘ E>1C i, TC iVQC'l 
CO MM on / lO VE R/ I OV ER 

Cb'MMON/XTJl/XTJl ■ ■ 

co hm on /F en CE /I fe nc e. AF en ce .b FE NC E. xF en ce 

qi me ns ION ' WT MOLE' (i O)' ' - 

D I MENS ION • ERAD (2 0) 

dimension TYPE2(¥)^' 

dimension TYPE3(2) 

DI mens! ON TYPE 9('2) i PHI PRT2T '' 

data PhIPR'/3H ,3HPhI/ _ „ 

da TA" T YP E9 /3 HC PX , 3 HP H I /' " ' ' “ " 

data ' TYPE3/3H W_,3HGAM/_ 

' daY7‘ tY'pe'? /iYh' '^ ozzlY' YT ohcent 

data XXP/l.E^-06/ 

■" "■■DATA'''iSHOC7-10 00 0/'' " 

WRIT£(6.40Q) 

■40 0~FO'RM AT (IHl ) ■ ■ “ 

j= 0 TWO oi mens IONAL 

J5Y' AX IS'YM METRIC''” " ■ 

SPECIES 1 IS H, . 

‘SPEC lES ■ 2 IS (5‘ 

SPECIES 3 IS H20 . _ _ 

SPECIES 4 IS" H2 ■ " ■ 

SPECIES 5 JS 0_2_ 

SPEC Iirs'''6' IS'"OH 

SPECIES 7 IS N2 _ 

species' '8 IS C02 

: rAV-..- 


c SPEC lES 9 IS CO 

c SPECIES lb IS C3HQ 

lPT^=q_ 

NX XJ =2 

IFLIP=0 

ihalt=o 

ICMPLT=0 . 

r" nstar=o 

I0PUT=32 " 

IEXIT=Q 

1AA=1 

"call I no at a 

XCOWlH=XCOWL 

' IEQ5=IEQ+1 

YN0Z=Y8P ; 

PEN= * PIN 

XJlSv=XJl _ _ . 

Pf IN F= PI NF «1 44 ;/ PI N " 

1FENC1 = 0 

■7211‘ IPP=0 " ■ r ■ “ 

iSTTrO 

WR IT”Er6» WO) 

IF (I FE NC 1. EQ ,1 ) WR IT E( 6i 31 00 ) _ ' _ _ 

"310 O' FORM at (3 OX ENCE INITIAL PROFILE*) 

IF(IfENC1.eQ,1) go To 31Q1 
“IF (I PTP. EQ .0 ) WRI TE (6 .9191) XOR(l) 

9191 format ( 2 ox«i Ni ri al Profile* 

-'"OR IG IN'* of SYSTEM =*E 13“.5 ) ' ' ' ” 

1F(IFLIP,EQ,0.AND, !PTP,EQ,1. aND.XFjD.NE.XCOWLH) write (6,919 2)' 

1 X0R{NXXJ-1) 

9192 FORM AT ( 16X«CHANGE OF ORIGIN PROFILE* _ _ ^ 

1<^ - ORIGIN OF SYSTEM =*£13.5) 

IF (I FL IP ,EQ. 0. and. IPTP ,Eq. 1. and. XF (1 ). EQ .XCOWLH) W R I TE (6 ,3 13 2) 

1 XOR(NXXJ-l) '*• “ 

3 13 2 FORM at ( 3X*PR0FIlE F OR U nD ER £X PA ND ED 5 HO CK _A T_ CO WL *_ 

!*■ '• ORIGIN OF system =*E13.5) ' ' " ' 

IF (I FL IP .EO. 1) WRITE(6,6oO) XOR(NXXJ-l) 

60 0 FORM at ( 8X*PR0FIL6 FOR 0 VE RE XP AN DE D" SH OC K" AT COW'L* 

1* • ORIGIN OF system =*Ei3,5) 

■‘rrn'FL IP .E'Q. 2)" WRI TE (6 ,621 ) XOR('‘NX'XJ-i') " " ” 

621 FORMAT ( 6X*PR0FILE WHEN SHOCK FL EC TS _0_FF LOWER_WALJ_* 

1 4 -" --ORIGIN OF SYSTEM =*E13.5) 

3101 CONTINUE 

IEQ1=IEQ-^1'' 

WRITE( 6,6896 ) PHIPR( lEOl) 

6 89 6 FORM AL ' ' ' ' / /5 X7*P T . * , 6 

IX, *X*, lix, *Y*, lix, lix, *T*, lix, *P*, lOX, *TH*10XA3) 
IFdCHEM.EQ.O.AND. ItQ.EQ.O) WR I T E ( 6 , 77 01 ) ‘ 

7 701 format (13X*ALP (1 ) = H*4X*Al.P (2 ) = 0*3X*ALP (3 ) = H20*2X*ALP (4 ) = H2* 

13X*ALP(5) = 02*3X«ALP(6)=DH*3X*ALP(7 ) = N2*) 

KSIDE=0 

■ INDl 3=0“ ■■■ ■ ““ ■ 

__316 0 CONTINUE 

FX2=0. ■ 

XLY2=0. 

‘ X MOM 2 = 0"."“ “ ■ 

XM AS 22 = 0 . 

*■ bo ■6 89 7 T=r»NP'TS ' ' 

YFPR1=YFPR 

THFPRI^THFPR — 
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YFPR=Yr( I ) 

thfprcthf( 1 )'” ■■■■ ■■ ■■■ 

IF (IfL ip .PQ. 1) YFPR=-YF( I) 

ir { 1 FL ip".eq.‘Y) thf'pr=-thf(T) 

_ IF (X Jisv .Eq. 0 . ) 60 To 3102 

IF (IPTP. EO .or GO" TO 3103 

.. IF (Xf( D ,E0. XCoWlH ) GO Tq 31 03 

IF(IFENCI.EO.I) GO TO 3103 

GO To 3102 ^ 

310 3 DO 3ill"j=i7NSP 

3 11.1,„.ALPDuM {J )= aL PF (J ,I ) - _ _ „ 

CALL ALL (CPXF { I ) Vw"F(T) ,R HO F'nT /RF ( f)TGAMF( I )"iEMF Y’l )7YMUFTI )7 

1 OF ( I ) , TF ( I ) , PF ( I ) , AL PD UM > 

IF (I .EO. iT GO" TO' 3”lb'2 ^ 

JQ =1 -1 

DY i3"=Y FP R- YF PR I " 

DXl3=XF( I) -XF( JO )_ _ 

THS13=1V57 07 96 3 

_ _ IF (DX13. 6T .1 .E-IG) T HS 13 =A TA N( DY 13 /D XI 3) 

2i=AXX (NXX j-iHK XF (JQ) -XBP )+8XX( NXXJ-1 r ■ 

Z2=AXX (NXX J-1) XF (I ) -XBP )+8XX( NX X J -1 ) 

XB AR 13 =■ ( XF ( I )' + XF’ ( J'O ) )'/ 2'.' -X SR F t " ' 

YB ARl3 -( YF ( I ) + YF (jQ) )/2. -YSHFT 

E1=G AMFCJQ )*PF (JQ) <^EMF (JQ) »EHF UO) <jSTN (THSr3-THFPR 1) /STN (THS13') 
E2 = GAMF(I )«PF(I ) <^EMF ( I _) «EmF ( 1 ) I N ( THS13-TH[f^R )/SIN(ThS13) 

F1 = (E1<*C0S(THFPR1)+'(PF( JQ)-PlNy)'' ' 

F2=(E2i>C0S( THFPR ) + (PF(I )-PlN)) 

El 1 = R H 0F'{ J Q r« Q f"( J o') « S I N (Th 5 13 - T H”f P R 1 ) /'S I N ( T H S 1 3 ) 

E22 = RH0F(1 )«QF( I )«SIN(THS13-THFPR )/SlN(THS13) 

ZBAR = (Zl+Z2)/2.' ' 

XMS13= (E11+E22 )^»DYl3/2 .-^ZBAR 
"XHAS22 =XMAS2?+XMS13 
FF X2 = ( FI +F 2 > /? . * DY 13 >»P F I NF »Z B A R 
. CO TT'HS =C'OSirHS 13')/ ST NrTH“Sl 3) ‘ 

XL LY 1= (E 1* Si N( TH fP R1 )+ (P f( jO )- pl N) QT TH $) 

XLLY2= (E2^S1 N( THfPR )+(PF'(i )• PFN) OT'TK S) ' ' " '' 

XXLY2= (XLLY1+XLLY2 )/2. <^DY13»PF INF 8A R 

XHOM 2= XhOM 2-f FFX2<fYBAR15- XXLY 2^^XB AR13 ‘ ‘ " 

FX2=FX2+FFX2 

XL'Y2'=X LY 2+ XX LY'2‘ ^ 

IF(I .LT.NPTS) GO TO 3102 _ 

IF(IPTP.NE.O')'"GD'TO 3105 “ ' 

FX1=FX2 

■“ 'XLY1=)<LY2 ■ ■ ‘ ■ 

XM0Mi=XM0M2 

XMAS"12 = XMAS22” 

GO TO 3102 

3105 TF(lFENCl.EQ.l) GO T0~3l06 ~ ' ' " ’ ' 

IF ( XF ( I ) .NE , XCOWLH ) GO TQ 3102 

3106" ■XTHXMS=FX2-FXiLxTHX ■ " ‘ " " - - - 

XM OM MS ^X MO H2 -X MO Ml -X MO M 

RA TM =YMAS2 2/WAS 12" -- - - 

EMS=RATM-1. 

IF ( ABS (EMS ) . LT il .E"-03)''“60" T0 ''316i ^ . 

IND13” IND13+1 _ _ _ _ 

IF ( IND13.GT. 1) GO TO 3l6l" ‘ 

CALL 'XMASSS(RATM,NPTS) 

I - - - - - _ 

3161 CONT INUE _ ; 

kside=i ' ■■’■■■■ ■ 
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3102 CONTINUE 
6897 CONTINUE 

DO 3i6 3 1 = 1# NPTS 
C=PF(I)7PiN 

YF PR = Y F( I ) 

■ ■ THFPR=THF( I) ' ’ 

IF (I PL IP .EQ. D YfPR=-YF( I) 

; IF(IFLIP.EQ.I) THFPR=-THF( I) ' ' 

IF(IeO,EQ.1) go to 9806 . _ _ 

WR lTE( 6r 16 ) 'I, XF (I r ,YFPR \ Qf (IV^fF ( I ) , C* 

WRI TE (6 ,160 2) (aLPF ( J, 1 ) , J =1 ,N SP ) 

GO TO 3163 ■ ■ 

9806 WRITE(6i16) I , XF ( I ) ,YFPR , OF (I),TF (T),C ,THFPR,RF(I) 

316 3 continue ” " "" — 

1800 format (« UNDERSURFACE AND CQNL THRUST =» El 3 , 5» 7X «L IF T =^^E1 3 . 5, 7X «M 

lOMENT =<^'E13.'5) . ' “ ' — 

__ IPTP=1 __ _ _ 

‘xjisv=xji ’ ■” ^ 

WR ITE( 6, 3130 ) 

313 0 FORMAT!/) “ - 

IF (KSI DE .EG, 0) GO TO 3104 

WR IT E( 6T18 00 ) 'XT HX ' .yLF T, XM Oh ^ - 

WR IT E( 6. 31 07 ) XT MX MS ,X hO MM S 

3 10 7 format (5 X^ SI DE WALL THRUST =« El 2. 4, 5X ID EW AL L MO HE NT =«E12",4/7) 

3 lO 4 CO NT iN UE 

LMAXrl 

" KMAX=2»(N”1) -I"' 

IF ( I OVER .NE. 2) 60 TO 7423^ 

IPP = 1 

LMAX=2 ■■■■■ ■ ' * * * - . 

km AX =2 

7 42 3“ NN1=N-1’ . 

DO 7 42 4 1= If NN 1 

■ II =NNi-i +1 ~ : 

X( I) =XF (I I) 

Y( I) =YF U 1) ■ ■ 

^ Q{ l)_=OF (I I) 

■ T ( I )"^T F* ( I I ) ^ 

P( I) =PF (II) 

TH(I) = THF (II) — — 

W( I) =WF( II ) 

R( I) =RF( II ) ‘ ■ 

DO 505 J = 1,NSP 

505 ALP( J, I)=ALPF (j/lT)'' “ 

DO 1 50 0 J=1,NSP 

150 0 ALPDUM ( J ) = ALP( J, I ) ' ~ ' ' “ - - 

call All (Cpx ( I > , w( i) ,rho ( i >, r( i ) ,g am ( i >» em ( i >, xmu( i ) ,q ( i ) » t( n , 

1P( I) , ALPDUM) - - - - 

7424 CONTINUE 

50 0 CONT INUE 

IF (KmA X. GP .1 91 ) GO T 0 91 23 

6060 IPP=IPP+1 ■ ■ ■■ - . -- .. 

6883 wR ITe( 6» 7633 ) iPp 

7 63 3 FO RM AT ( 5 X , «-L I N E NO ^ 1 4 ) ’ * ~ “ ' 

IF (I OVER .EG, 2. AND. KMAX .GE. 30 ) GO TO 67 0 

JJ=i+KMAX/lO • - — “ •- 

JJ1=KMAX-1 

IF(IO.EQ.I) JJ = 1. - - - 
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__ WRITE(6,6885)TYPE3(IEQ5) »TY_Pe9(IEQ5) 

-16 rORf-1AT(3X.I4,3X/lOE12.4) 

6685 rORH aT (5X» «PT. «_i 8_X ,HX »_^Y_« »llX , . 1 pX > ^ TH > 11 X > 

' '‘llOXi‘«E'M'^ iroX7A‘3 V 9X/A3) 

DO 7 63 7 I=1»JJ1^JJ 

I)/P1N' ' 

YP pR =Y ( I ) _ _ _ _ 

THFPR=TH (ir 

DUpR =W ( 1 ) 'V 

I fT\ EQ 5 . EQ , 2 ) DU PR“=G ‘aM TI ) 

cPP=cPx(0 ■■ — 

“ IF ( lEQ .EQ: i) CPR-:=R ( I ) 

IF (I fL IP .EQ^ 1) YpPR=-Y (!) 

IF“(IFL'IP”E0. irrHFPR=-TH (I) 

WR 1 T E ( 6 » 76 32 ) U X ( D , YjP ^HF jR^QjJjj^T (JJ_JJlU 

7 63 2 FO RM AT (3'X» T473 X, lO'El 2. 4T 

7 63 7 cO nT in liE _ : — 

1 60 2 FORMAT (1 OX, lOE 12, 4) 

jF { I pP .E Q. 1) _G 0 TO _7 75 9_ 

"r=KM AX' ■ 

C=p( D /P IN 

YF PR'=Y ~(T)” “ 

THFPR=TH J IJ 

DU PR ■( I ) “ 

IF(IEQ5.E0.2) DUP Rj G A M_( 1 )_ 

CPR = CPX'n) 

IF(I£O.EQ,l) CPR=R(I) 

I F rfrCl P^ E Q , 1)" Y F P R = - Y ( I ) 

IFdFLIP.EQ,!) THFPR=_-TH J U 

"WR ITE ( 6, 76 32 ) 1,X(I),YFPr,C , t HF PR , Q ( I ) , T ( I ) , EM < I ) , DUP R , CP R 

___ GO TO 7759 

670 'JJ=1+KMAX/10 

"“-fs'H= IS HO"C ^ ™ 

JJ1=ISH-'1 . 

WR ITE( 6, 68 65 )T YP E3’(1 £Q'5T 

67 2 DO 6 71 1 -<J J2 , J Jl , J J — ~ — — 

C=P( 1) /P IN “ ^ 

YFPR=Y ( I ) , ^ 

THFPR=TH"nr 

DUPR = WCI) _ __ ___ 

1F ( IEQ5.ED";'2)“DUPR=GAM('nr 

CPR=CPX(D _______ _ -- — 

IF(lEQ.EQ.irCPR=R(I ) 

lFjiFLlP_^E_0.1)__YFPR=-Y (H 

TF fl FL’ IP . E 0 .'l ) T HF PR = “ t H ( I ) 

WR ITE( 6, 7632 ) I , X( I) ,YFPR,_C ^ » T HF PR ,_Q77 K T^^ 

-■’671 CONTINUE ‘ 

I = JJl+l 

C=P< I )'/P IN ■“ 

YFPR=Y ( I ) 

THFPRVtH'Tir 

DUPR=W(I) _ 

■■ IF(IEQ5.GQ.2r'DUPR=GAM( i ) 

CPR=CPX( I) ■ . 

■ IF ('I EQ .EQ, 1) C;PR=R ( I') 

IF( IFU IP^EQ.l) YF_PR = -Y_<n 

jFdFLTP.EQ.ir THFPR = -tH (I) 

WR ITE( 6, 76 32 ) I, X( D ,YFPR, C ,THFPR *0. f d aFMP_vOUF_R' OPR 
IF ( JJi .EQ. KMAX-1) GO TO‘ 7759 '' 
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JJ1=KMAX-1 

GO To 67 2 : 

7 75 9 BE TP R= BE T 

WR ITe( 6> l8o0 ) XThX ,YlFT. XMOM _ _ 

“■ IF ( I FL IP .EQ, 1) BETPR=-BET 

IF (I over .EO. 2) WpITE(6 ,653 ) BETPR^EMCI 

653''F0RMAT {5X>-^SHOCK ANGLE =«E12. 4, 5X^^EXTERNAL HACH NO. r*E12.4) 

WR I T E ( 6 » ^6 03 ) 

160 3 FORM AT (/// ) 

6684 iCMP=IcMPLT+1 _ 

IF ( I OVER .NE. 2) GO TO 2656 

ISH0C= lSHOC+1 

IF (NST AR ,EQ, 1) I SH OC = I SHOC-2 

2658 CONTINUE ' 

■■ IF'( fHA LT‘;E QViT ’S Td P' " 
if(1flip.eq.2.and.n.gt.npts) stop 

GO TO C4 30 0, 43 01 » 4 30 2, 43 03 )‘, ICMP 

430 1 IF < I CO WL .EQ. 0) STOP ___ 

ICMPLT=2 ' ■ ■ . ■ ■" ‘ ' ■ 

CALL C0WL(-1.) 

'iPM=r 

4 30 2 IPM=IPM+1 

IF (IPM .GT. MM-3 > GO TO 43 03 ' ” 

LM AX =k ka x^- 1 

X N (1 )= X M (I PM ) ‘ ■ ' 

Y N (t )= Y M (I pM ) 

p - N (1 ) = -p- MM (] PM ) — 

W N (1 ) = W M (I PM ) 

0 N(1) = Q M(IPM) ' “* 

R N<1)=R M(IPM) 

T N(1)=T M(IPM) ’ ' ' " . ' 

IN N(1 )=TH M(IPM) 

■■■ ■■ EM N (l')= EM M (I PM ) ' ■ 

RH ON (1 )= RH M (1 PM ) 

GAMN a )=G M (I PM ) 

_ XMUN (1 )= XMUH (I PM ) __ _ _ * _ 1 

CPXN(1)=0. ■■ -- - - - 

IF(IEO.EO.O) CPXN( 1) =CPXM( IPm) 

“DO 4 304 J=l. NSP 

ALPN(J,1)=0, _ 

4304 IF ( 1 EQ .80. 0) ALPN( J» 1) =ALPM( J, IPM).' 

GO To 80 60 _ _ _ _ 

^303 ICMPLT=3 

STOP _ _ _ __ _ 

4300 IF(X(i>.EQ.XCOWLH/AND.TPP;EQ.NPTS'.AND.T6VER\NEi‘2) GO" TO 4301 
IF (N .GT. NPTS ) GO TO 8104 

“ LMAX=LMAX + 2 “ ■ ' ■ “ 

L=1 _ _ _ _ _ 

K=i ■■■ ' ■ ' ■■ 

_ DO 510 J=1,NSP 
'5i0 ALPN(J ,L ) = ALPF ( J,N) " 

XN<L)=XF (N) 

YN (L ) = YF ( N) ' 

QN(L)=QF (N)_ 

' PN(L ) = PF“'(N) 

TN (L ) = TF ( N) _ 

■ THN( L) =THF '( nT 

WN(L) = WF(N) 

RM(L) = RF(N) ■ ■ 
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DO 1501 J=1#NSP _ „ 

i50r ALPDUM ( J ) = ALPNi: J', 1 ) 

CALL ALL ( CPXN (1 ) , WN U ) , RHON (i ) r RN (1) » GAMN (1 ) # 6MN( 1) ,XMUN( 1) i QN( 1) 

tvtn( i'r»pNri)';ALPDuM)” 

_ IFeNc1=1 _ _ _ 

ir'(rFENCETEQ‘:rrA'NOTXN'( LT.EQ/XCWLH GO tO 26 01 

irENcl=0 _ __ 

'“GO'TO 8060 ■ ■■■■ ■ 

810<5 ASL = TAN(TH(1)) 

■" A C H ='T A N ( Th“( '2 ) '+ X M U ( 2 > ) 

_XN a )= (Y (1 )-Y( 2) +XJ2 ACH“J(( 1) <»ASL )/_(_ACHj-_ASL )_ •_ _ 

DELTH=XN(1)“X(1) * 

_ IF (( DELTH/TEST )_^GT .1 .J GO TO_ 9 06 1_ 

ICK= 1 

_151 Hf_Y ( 2) +A CH » f_X(^( !)• XiZ) ) 

C A LL '8 OD 0 ( XNl 1 T, YN (1 ") . TR'tMT 1 ) ) 

__ IF (I OV ER ,Ne, 2) Gq_T0_6_94 

YN (1 )= -Y N( 1) 

THN( 1) =-TH_N(_l) _ 

"69 4 CONriN UE' ' ” T 

DELTH= ABS( THN( 1) -TH( 1) ) 

■ IF (DFL‘TH7GTV ro'2)"W*‘TO "906"r 

ER=( YCH“ YM (1 )) _ ‘ 

'iFCA'e'S (ER) ,LT/1, H-0'4 ) G'O TO 3510 

^icK=icK+i ^ : 

' IF( ICK.LT.15) GO TO 3511 

WRIT£(6,3512) _ 

■35r2“''F0RMXTTi>“"Td0l/iANY 'I N Tn I f I AL’"' BO'DU CALL FROM NOZBOD<») 

STOP _ __ _ 

”35il ’ IF ( TcK ,Gf. 2) GO TO 3513 

ER1 = 6R 

XSL1=XN(i) 

XN (1 )- 1. 01 «XN( 1) ; ^ 

. GO'TO' 35 K" 

3 5l 3 rii H= XS Ll “I Rl XN (1 )- xS Ll )/jE R- £R 1) _ 

■ ER1=£R "■ ■ "■ 

XSLl=XN(l) 

"XNU')=DUM ■ " 

GO TO 3514 

T51T C0NTTNU£ " "" - - ^ 

XSL=XN(1) 

VSC=YN(1) ‘ 

THSL=ThN(1) __ _ __ 

lP(roELT'H/TEST ) /CTVi:T GO T0“9 060 - -- 

_DELTH=XN (D^XU) 

9061 LMAX=KMAX+1 

DO 9 00 0 M1=2,KMAX_ 

L=KMAX-M'1+2 ' 

i^*"i _ _ : 

X ('D =x" a') ' 

Y (I)=Y (L) 

I)'=Q"" (L) '■ '■ ■■ " 

p ( I) =P (L ) __ 

' T (D=T’ (L)^ " "" . ■■ ■■ '■■■" 

W (I)=W (L-) . _ __ _ _ _ __ _ _ 

R ( I ) =R ■ (L ) 

TH_ < I ) JL ) ^ 

EM 6r)‘=EM’ (Ll 

CPX(I)=CPX(L) _ ____ _ _____ 

RHO( I ) =RHO (L > 



GAH( D =G am (L ) 

XMU( I) =XHU (L ) 

DO 9 OO 1 

9 001 ALP( J» I) =ALP (J #L ) 

_9 oO 0 cOnTiNuH _ _ 

ISHOC= ISHOC+1 

kmax=kmax^*i _ _ 

, RAT=.5 

__X__(2)=X (1)+RAT«(X (3)-X;‘(i)) 

Y (2)=Y (i) + RAT«(Y (3)-V“(l)) 

0 (2) = Q (1)+RAT<^(Q (3)~Q (D) 

T (2)=T (1) + RAT<^(T (3)-T (D) 

TH(2) = TH(l) + RAT«(TH(3)-THa)) 

R (2)=R (l) + RATi^(R (3)-R (1)) 

CALL P AN DW ( P ( 1 ) ,_P ( 3 ) _^R AT 2 ) , W ( 1 ) » W ( 3 ) » Q Q(3)iO(2),W(2), IEQ»0) 

DO "9 00 3' J=T, NSP ■ ~ ' ' ■ 

ALP( J» 2) =ALP (J ,1 )+RAT<^ (ALP (J ,3 )- ALP( J» 1) ) 

'9 00 3 ALPDUM (J )= ALP( J, ?) ' 

_ CALL ALL (CPX (2 ) , W( 2) fRH0j2 ) , R{ 2) ,G AM (2 ) . EMJ2 ) , XMU( 2) ,0 (2 )» T( 2) , 

“ '1P( 2) , ALPDUM) ■ ■ ' ■ 

GO TO 8104 

”'9'06 O 'LM AX =K'MA X ” * - - — . . 

A"l. _ 

B=0. ■ - - - . -- - -- - 

IF (I SHOC .EO. 1) 8ETN==BET 

■’8 00 0 CO NT IN UE ----- -- 

ThN(1)=THSL 

8 03 0 EMSL =XM3 ( . 5» .5 ,TH( 1) ,ThN (1 )) 

EM1 = XMK Ai B» TH(2 ) » XMU(2) ,THN(i ) , XMUN(l) ) ~ ' ' 

IF( ISHOC.EQ.I) EM1= .5»(TAN(BET)+TAN(BETN) ) 

XN(l)=(Y(l)-Y(2)+X(2)<»EMl-X(l)»EMSL)/(£M3«-eM$L) 

__ICK=1. 

3 52 4 YC Y( 2T +E Ml XN (1 )- X(‘2) > 

_ CA LL B OD U( XN (1 )» YN (1 )# TH n( 1) ) 

IF (I OVER .N£. 2) GO TO 6 98 " ' ■"“ '7 

YN (1 )= -Y N( 1) 

THN( 1) =- THN( 1) “ 

69 8 CO NT INUE 

ER=tYC H-’YNTin ■ ~ 

IF(ABS(ER).LT,l.E-04) GO TO 3520 

ICK=ICK+1 ““ - 

IF (I CK ,LT. 15 ) GO TO 3521 

WR ITE( 6/3522 ) " ' 

_3522 F0RMAT(« TOO MANY ITERATIONS IN FINAL BODU CALL FROM NOZBOD») 

STOP ■ '■ 

3521 IF (I CK .GT. 2) GO TO 3 52 3 

ER1=ER ■ " ■ ■ ■ 

XSL1=XN(1) 

XN (1 )=1* 01»XN( 1) " - - 

GO T 0 35 24 

3523 DUM=XSL1-ER1«( XN a')-XS'Li )7TER-ERD" 

ERl=ER 

XSL1=XN( 1) ■ ‘ ■ - - -- 

XN { 1 ) = DU M __ 

GO TO 35 24" - 

_3520 CONT INUE 

"8020 7i=f k a, b7xmu( 2)Vg am'"(‘2 )TpT2T7xmun'( 1) ;g‘amn( 1) /PNTi ) )‘ 

XDUM2=X( 2)“XCR (NXX J-1) 

X.r)U.'lN=XN (1 )-XOR( NXXJ -1 ) ' ' ~ ' “ • 
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A2=F2( a. B. J_, Xvilj XDUH2»J( 2) »Th( 2) |XMU (2 U XDUMN_, YN' a 

'iN(i)') ■ ■ ■“ 

__ DUMl =A 2 <mXN( 1) -X (2 )) 

1F'(X jf.GTro.') ”DUM1=A2^AI: OG (XDU MN/XDUM2 ) 

_ __DUM= _ < TH (2 )- THN( 1^-DUH) /A 1 

PN (lT= PT 2y «EXPlnUM ) 

IF{B.E0.0.) GO TO 18 01 _ „ 

■■■■“ 

^L2=l 

TpYX isl(Tl ) . GT . X CO WL -1 . E '"D 4 ) GO TO 1801 

IF (X N( LI ) . GT .X XJK NX-XJ )- 1. fe- 04 ) GO TO 18 01 

Zl=AXX (NXX J-1) »( X( L2')-'XBP) +BXX (NXX J-1) 

_ Z2=AxX (NxX J~l) XN (LD -XbP )+BXX(_NXXJ-l_) __ 

DAX= (Z i+ Z2 )/2. «A BS (YN( Ll )-Y( L2 ) >“ 

da y=_(Zjl+ z2 )/ 2. IL^JJ 

■>0'r (PN( LI )+P( L2 ) )'/2 . 

„ PA V= (P aV /P JN -1 , ) IN F*_14_4j_ „ _ 

DX TH X= PA DA X ' 

dylft=-p AV^jDaY _ , 

XHS= (XN( LI )+ X( L2 )) /2 . - XS HF T 

YMS= (YN( Ll )+Y( L2 )) /2 YSh_F_T 

'dmoh =y>;s XT Hx -xks«d yl ft 

XTHX=XTHX+DXTHX 

- Ylft‘=YLF t+ dylf t 

XMOM = XMOM+nMOM . 

iWl ‘ CONfl NU’E " 

DO 8 0 5 0 J^ l N S P^ 

"a irPNTJ » 1 ') = A L P 7 J 1 ) 

e_050 ALPDUM(J) = ALPN(Jrl) „ 

IFUEOTEQ,!) go to 3000 

TERH2 = RH0(l)i^Q(l)<^A _ _ 

‘"'if (B .GTTO. ) TERM2= TERM2+Bi^RH0N ('l j>tQN (1 ) 

IF(B.E0, 0.) CPXNji )=CPX( lY 

. QN''(r')= '<P (i )- HV (i )'^ /T ER M2 +Q a, ) 

TN (1 )= T( 1) +( QN (1 ) + Q( 1) )* (P N( 1) ~P (1 )) <»E IN / ( CP X( 1) +C PX N( 1) )/ TE RM 2 _ 

qA-LL~AlL (CPXN( 1)VwN{ 1) ,RhON( 1) /R N( 1>'7G AKn(I) ,E HN (1')VXMUN'(1 )i 

IQN (1 ), TN (1 ), PN_(l ), ALPDUMJ ; 

" GO TO' 30 01 

3 00 0 CO NT IN UE . __ _ _ 

RHON'a )='RB O( r u ('P'Nrrr/PTi ‘i ^vir/G'AH n > > 

VVL=Q(1 )»Q(1 ) + 2. -»GAM a^)/(GA_M( 1 )-l.)»(P(l )/RHO(l )-PN(l)/RH 

iON(l))" '■ ' --- -- - 

qN(1)=sqrt (Vvl) _ _ ; 

WN U' ) = W ( 1 ) +'. 5 «■■( Q ( i ) « 0 (”1 ) - V V LT' “ 

jn(1)_=R(1) „ : 

"" T N i’ ’) = FT ( P N ( 1 ■) 7R N ( 1T , W N (1 ) ■) 

GAMNd ):;FGAH(TN( 1) .PN( 1 ) ,RN( 1) ) 

* £MN(1) =ON( l)/SORT(GAHN(l )»PN(1 )/RH0NU-)') 

xmuN(1) = ataN(1./sqrt (Emn(1) ««2'1 . . : ^ 

3001 CONT INUE 

1 / ^ 

a=:5 

0 = ^5 

' ‘ GO ‘TO ' 80 00 ” ^ " 

6 049 IF ( I SHOC .NE. 1) GO TO 2600 

call SH0CPT( ISH0C/3»8ET, BE1N,A ) 

BE Y= TH N (_2 )_- (J3E TN -TH N(1 ) ) 

2625 CALL^SH0CK(BET,QN(2) rTHN(2)»GAHN(2),EMN(2) iRH 0N(2) ,PN(2) ,NN(2) r 
IRN (2 ) , QN { 1 ) » TH2» GAMN (1 ) ,'EMN ( 1 ) » R HO N ( 1 } # P N ( 1 ) , WN( 1) rRN( 1 ) dN( 1 ) , 



1 (1 ) » -1 . > _ _ 

ER4=THN( D-TH2 

(ArS (Er4j^ LT ,1 .Ej-Oj)_Go_TO 2 60 l 

lGG=lGG+i * 

IF(IgG.GT.15) go to 2627 _ 

IF(1GG.GT.2) GO TO 2628 ' ’ ‘ 

_ ER1=ER4 __ _ . 

BETP=8ET ■■ ■ '■■■ ‘ 

BET=8ET^ ,02^ ^ ■ 

GO TO 26 23 ' 

2 62 7 WR ITe( 6» 26 29 ) 

2629 format ERROR IN BETA LOOP’lN M’ai N«T‘ 

STOP 

2626 DUMM=8ETP-ER1« (BET-DETP) /( ER4-ER1) 

ER1=ER 4 
BETP=9Ef 

eET=DUMM 

GO TO 26 25 - - 

2 60 0 CCMT INUE 

IF (I FL IP ,E O/i; AND;XN (i ). GT .XFi NAL) IHALT=Y ' 

JF (XN( 1) ,LE r XCO'i^L- 1> E" 04 . A ND , X N ( 1 ) . L E . XX J1 ( N XX J ) -1 . E -0 4 ) GOTO 80 60 

XNN=XCOWL ' ■ ‘ 

IF (X N( 1) .G T. XX J1 (N XX J) -1 ,E -0 4) X NN =X XJ'K NX XJ ) 

“RA T= (X NN -X a )) /( XN (1 )~ X( 1) ) ' “ 

XN (i ) = XNn 

CALL 80DU( XN (1 ), YN (1 ), THN{ 1) ) - ' ---- 

p at m p Njij- p ( 1) ) _ 

L2=l 

Z1 = AXX(NXXJ-1)«(X(L2)-XBP)>BXX(NXXJ-1) ' ^ 

22=AXX(NXXJ-l)»(XN(Ll)-XBP)+8XX(NXXJ-l) 

DAX=(21+Z2)/2,<^ABS(YN(L1)-*Y(L2>) ' ■ 

DAY= (Zl+22 )/2. «( XN (Ll) -X (L2) ) 

PA V= (P N( Lr )■+ P( L2 ) ) V2 . ^ 

PA V= (P AV /P IN -1 . ) «P IN r» 14 4. 

DX TH X= PA V» DA X ‘ ; - 

DYlF T= ~p AV*DaY 

XMS= (XN{L1) + X(L2))/2Y ‘-XShFT ^ 

YMS= (YN( LI )+ Y( L2 )) /2 YSHFT 

"DM OM zY MS »D XT HX ~X MS «D YLTT 

XTHX=XTHX+DXTHX 

■ ylft=ylft^dylft ■ --- ■ ---- 

xmom=xmom+dmom 

0 N(1) = Q (D+RAT«(0 N(i)-Q (1))" 

T (D+RaT«(T N(l)-y (D) 

RN (r) = R{ i’) +RAT*( RN (1 )-R( 1) ) 

WNT= WN (1 )-»-QN (1 )»QN (1 )/2. 

WT = W (1 ) + Q( 1)«Q (1 )/2, ■ “ . 

WNT1=WT+RAT<^ (WNT-WT) 

WN (1 ) = W< 1) -fR AT »( WN (1 )- W( 1 ) ) ’ 

iFdEQ.EQ.l) WN( 1) =WNT1-0N (1 )«ON (1 )/2, 

DO 3 519 J=1,NSP' "■ . 

ALPN (J ,1 )= aLP( j, j ) 

3519 ALPDUM (J ) = ALP( J, 1) - - — — 

CALL all (CPXN( 1) ,WN( 1) ,RhON( 1) ,RN( 1) ,GaMN( 1) ,EMN (1 ), XMHN (1 ) , i 

IQN (1 ) , TN (i ) , PN (1 ) , ALPDUM ) ■ ~ “ 

2601 continue 

IF ( X J1 . EQ. 0 . , AND . I OVER , EO. 0) GO TO 6060 


-AIO" 



IF (I rL IP .Eq. 1 ) Ix=2 

2655 IAA=IAA+1 
XF (I AA )=XN ( I X) 

>FlI AAy=YN'( l xr 
PF (I aA )= PN < I X) 

OF {rAA ) = QN (I X) 

rF (I aA )=RN (I X) _ _ _ 

wr (I AA')=WN'(rX) 

TF (J aA )=JN (I X) ' 

THF(irA)‘=THN (TX)" 

XMuF ( I aA )= XHuN ( I X) 

GAMr ( I aa ) = gamn ( I X) ■ ' 

£M F ( I aA ) = EM N ( I X) 

“'RHOF'd AA )=RHON axr' 

DO 3 47 0 J= 1» NS P _ 

3 '47 0 Tl'PF' (J n'AA")='AL'PN'<J , I XT 

IF (1 X. EQ .1 ) go T o _26 56 

IX =1 

_GO TO 2655 ^ 

2656 IF (I FL IP .EO'.'i) I AA =fAA "i 

lChPLl:^0 

NS TA R'= 0 

IF (XN< a> +» OOOl .GE, XT Jl^ XJ1=0,_ 

• iF U FL IP'.NE, 1 

IF(IFENCI.EO.I) IFENCE=0 

"“1F(NXXJ. GT .NXXJl ) NXXJ=NXXJl+i 

I R XN (_1J . EQ . XCq_WL )_ XCOWL = XFJ NA_L 

IF (TflI P G T . 0 ) GO" ■ T 0 ""2 6"0 2 

IFdCOWL. .EQ,0) GO TO 7593 _ : 

RP-PMM (MM) /PF ( IAA)ttP IN 

IF (RP . GT . 1 . )_ GO TO 5989 ^ 

7593 NPTS=IAA 

1 0 VE R= 0 

*IA'A=1 ' 

GO TO 72 U 
5 98 9 CONT INUE ■ 

EMEX T= fcN INF^^QM {N m) /SqR T( TM (M tA) ) __ 

" CA LL C OW LO ( I AA , T HM ( M M ) , E ME XT , B ET' , P'mM (”M M F) 

XI NTu=X^ - 

"DO “6 91 fl3'=l ;3" " 

XXU( I13)=XXL(I13) 

AAK I13/1T = Aai ( 113 72 ) " *' 

AA2( I13>1> =AA2 (1 13 ,2 ___ 

■691 AA3('I13R)=AA3(I13V2) 

2 60_^ cont inue 

NPTS = IAA+i 

IFL1P=IFLIP+1 _ _ 

IF ( I FL ip .GT: 1) ■XXP=XFI NAL ' 

IAA = 1 _ _ _ 

LSTT::0 

I0VER=2 

call rNDATA * 

BET = -BET _ __ 

ISHOC=l " 

GO TO 7211 ' _ _ 

8060 continue 

LM=LMAX-1 

I F ( X J1 . G TV O', UM ; I AA , G T . 1 ) L M^lOOO 

IF (I OVER .EQ. 2, AND. L5TT ,EO, ly LM=1000 
"'"LST = 0‘ 

-All- 


M=0 

iDpTrO 

XXP=XFINAL 
IF (NsTaR .EQ. 1) M=2 
IF(ICMPLT.E0.2) M= 1 
5520 L =L'*'1 

'lF(L.GT,LM> GO TO 6520 


67 8 M=M+1 

g- Q ^ 

250 CONT INUE _ _ 

IF(L ,E0. ISHOC+l) CALL SHOC PT ( I SHGC'^M .0 ET ,B £T N, A , B ) "'■ 

IF(L.NE. ISH0C«^1.0R.B.EQ,0. .OR, IFLIP.NE.1) GO TO 2603 

IPL=3 ■ ' 

2684 P0LYC=-Aa3 (I PL ,1)- YN (L )+ TAN( BETN (XN( L) -XINJU ) 

P0LYa=“AA2(lPL ,1)-TAN(BETN) ' 

POLY A=-AA1 ( IPL ,1 ) 

IF (POL YA .NE. 0. ) ’ ’ ' 

IXXP =( -POLYe-SQRT(PQLYB»POLYB-^. -^POLYA ^P OLYC ) ) / ( 2 . «P0L YA ) 

1F(P'0L YA,E 0,0. T XX p POLY C/ POLY B 

XXP =XXP +XINTU 

IF (X XP , GE .XXU (I PL ))■ GO TO 2 60 3 

ipl=iPl-i 

■ GO TO 26 84 . - 

2 60 3 CONT iNUc 

IF (L .60. ISHOC+1) GO"TO 3 003'~ - .. 

EM1= XMK A, B, TH (M ) , XM U ( M ) , T HN (L ) i XM UN ( L ) ) 

EM9 = EM1 - - * 

IF(L,EQ. ISHOC) EM1=TaN(BET) 

IF (L .60, ISHOC. AND.H. GT .3. )EM1 = TAN (8ET + THN (L«1) + XMUN (L-1 I-TH (M-lV 
I-XMU(M-I)) 

EM 2= XM 2( A", B/ TH N( KT/X'MU N( ‘Kr,T HN' (L ') /Xm'UN COT' ' 

XN (L )= (Y N( K) -Y (M eM 1» X( M) -6 M2 *X N( K) )/ (E Ml -E M2 ) 

■ YN (L )= Y( M) +E Ml »( XN (L )- X( M) ) " ; 

IF (8 .Gt, 0. ) GO To 6Bl 

IDROP=0 - - - - 

BETTzBET 

ir(M .GT. iTANDVrOVER'. EQ;2' ^ ETT=BE T+'TH N(X-X) MQ N'(X-* 1 ') -T H( M-l) ~XMU'(M"~' 
1-1 ) 

'IF (L .EQ. ISHOC-1) call DR TE ST (X N( L ) VY N{ L ) ,E Ml ,B ET T, ID RO P, +1 . M ) 

IF (L .EQ, ISHOC+2) CALL DR TE ST ( X N ( L ) , Y N ( L ) , E Ml , B ET N , I D RO P , -1 , M ) 

IF (I DROP .EQ. 0) GO TO 6 81 " - - - . - - - - . . . . 

LMAX=LMAX-1^ 

lm=l'm-i ■ ~ ^ 

IF (L .L T. IShOC) I SHQC =1 SHOC-1 

GO TO 678 ■ 

68 1 cONT iNuE 

IF (L .NE, ISHOC) GO TO 601 ' 

_ „ YSHOc= Y< M) 

“XSHOC-XI M) 

_ IF (L .G T. 2) Go TO 6 02 

■'YDUM=YN(L>-EM9«(XN(L)-X(M)) ' “ ■"* 

_ XDUM -X (M ) 

■RA T= ( Y DU M= Y ( M ) ) / ( Y NC 1 ) - Y ( M ) ) 

Y(M)=YDUM 

X( M) =XDU¥" "■ 

QKO=Q( M) 

Q (h)=Q (M) + RAT«(Q NXYT’-Q' (M))' ~ 

-AI2- 



_ __ T (M)=T (M) + RAT«^(T N( 1)-T <M)) 

R' (M)'=R (H)"+RAT<^(R '■ ■■■{Mn" 

CA LL P an DW (P (M ) , PN (1 ) , RA T, P( M) . W (M ) ,_WM (1 ). Q HQ> Q Njl ) > Q( M) (M )_> lEQ 



TH (m)=TH (M)+RAT<»(Th N( 1 )'• TH ( M)_) 

’ XM U ( M ) =X fiU (■ M 7+ R A ( X MU N ( l') - XM U Tm') ) 

GAM( M) =G AM (M ) + RA (G AMN^__1 )-GAM( M) ) 

bO'644 ’j = i,N5P 

_ 64 4 AL_P< J» M_) =A lP ( J ,H )+ ^ T^1?-LP Nj J» Jj VAL_P1 J 

"go T0''60’i ““ '■ 

60? CONTINUE ^ • 

EM 8= XM ?( .5 r. 5r”TH (M‘-i ) , XM lj( M- i) , T H ( M) r X MU {M ) ) 

YDUM=< EM9^Y( M) -EM8 «YN( L) +E M8 «EM9«( XNjL _)~X( £M9- EM 8) 

XDUM =X N( L) -( YN (L )- YDUM )/ EM 9 “ ‘ 

_RA T=_(Y DU M- Y< M) )/ (Y (M )j:^Y(Jl) ) 

YTM’) =YDUH"' 

x(m)=xqum 

''QHO=Q(M)"‘ 

Q ( M> =0 (M)+RAT<»(Q (M-I)-q ( M) ) 

f ■"■■( M) =T’ ” (M ) + RA Ttt (T " (M -1 )- T ( MrY' 

R ( (M )+ RA (P (M-I)-R _( Mj ) _ 

CA'LC' P ATJ DW (PlM ), P( M-i) »RAT ,P”(H“)'; W<’H> ), QHO, Q( M-1) (M ), W( M) 

__ 1 , 1 EQ , 0 ) „ „ ^ _ 

■ ' TH (MXr'TH (M) + RAT»(TH (M-i)-TH (M)) 

xmu(m)=xmu(M) + rat«(xmu(m~i)-xmu(m) ) : 

G A M ( H y = G A M ( M ) + R A T ( G A M ( M - 1) - G A H ( M ) ) " 

_D0 604 Jj=_l»NSP 

60 4 A LP (J* M“) = aUp (J»M)+RAT«(ALP(J,M-i)-ALP(J»M)) 

601 CONTINUE __ 

“A1=F1( A.BVxmIjcM) ,GAM(M),P(M) rXMUN(L) ,GA”MN(L) »PN(L) ) 

81 =F 1 ( A , B » XMUN ( K ) , G A MN ( K ) , P N ( K ) » X MUN (_L)_, G AMN ( L ) . PNyL.).) 

XDUMK = XN (K )-X0R( NXXJ-1 ) 

XDUMm= X( M) -XCR (NXX J~l)^ ^ ^ 

' : XD'UM L= XN" ( L ") - XO' R < NX X J"-'l )" 

A2 =F 2( A. 1. ,X j, XJ V XD UMj^,_Yj M) .T H( ,X MU |M U XD UM L» YN_(L_) » TH N( L) »X W 

1N( L) ) ”■ ‘ 

B2 =F2( A» -1 XJ ,X Jl ,XdUmK ,Yn( K)_,J HN ), XM UN )^UML^YN (L ), TH_N<; L) 

IVxmu'nc dt * * ‘ " ’ ' 

DUM1=- (A2-^B2)«XN (L ) + A2'“»X (M ) + B2»XN ( K ) 

IFIX JT7G rVO.lD UMi= -{ A2 L"OGTXD UML/ XDUhi'K) +r 2« AL OG (X OU Mu7x DUNK ) ) 

DU M = ( A AL OG < P ( M ) ) +B AL OG ( P N ( KJ )_+ JH ( M ) ~ TH N( K ) +D UM 1 ) /J_A1_+_B 1 ) _____ 
PN(L”)='EXP(DUM) 

DUH=ALOG(P(M)/pN(L)>ttAl 

DUM2=XN( L) -X(M ) 

IF (X Jl.NE. 0. ) DUM2=AL0G( XDUHUXDUMM) 

^yHN^'C) =dUM + TH ( M ) -*A'2»D0M2 ‘ 

. 253 EM3A=XH3(.5, ,5tTH(M) ,THN(D) __ 

- -eM3B=XM3( .5 i’.5,THN{K ), THNCL)') ' 

8M3= ,5>“ ( EM3A+EM3B) _ _ 

' TESTX= ABS( XN (K )-X(H) ) ~ 

test y= ABS< YN (K )- Y( M) ) 

IF (T EStX .GE, i. e-io‘) 

ISL 0PE= (VN( K) -Y (H ) ) /( XN (K )- X( M) ) _ __ __ __ 

;■ IF (TESTX";LT7i:E- 10) ‘SLOPE^l.EiO 

[ )(D = (YN(L)-Y(M)+Sl0PE<^X(M)-EH3«xN(L))/(SL0PE“EM3) _ 

IF (NST AR .EQ, 0 ) GO TO 6881 

. XX=A3S(1 ,"X(M)_^XN(K) ) ' 

-'jF'(Xx/LT';TobT) XD= (xnTk) +x cm ) ) /2 , 

6681 YD=^YN(L)+EM3<i(XD-XN(U ) • _ __ „ 

■ ■ IF (TESTX .LE. TESTY) GO TO 2250 



kIP=0 _ 

2230 RATD = ( XD-X (M )) /( XN (K )-X( M) > ' 

THD=TH (M ) + rAtD«( THN( K)_-T H( M_) ) 

EM3=XM 3r.Y/, 5/THD, 

XDT=XD 

XD = ( YN (L )-Y( M) +SLOPE*X (M )-EM3ttXN (L > ) / ( SL OPE-EM 3) 

YD=Y^'( L) +Em3«( xD-XN( L) ) 

IF (ABS < ( XDT-XD )/TESTX) .LT, .0 01 ) "GO TO 22 31 ' ~ 

KlP=KIP+i 

""if (K IP ,GT. 20 ) CALL ERROR (2230) 

GO To 22 30 

2 23 1 RATD=( Xn~X (M )) /( XN (K )-X( M) ) ' 

GO T 0 22 35 _ 

2 25 0 CO NT INUE 
KIP=0 

23 d RA TD =( YD -Y (M ) ) 7( YN (K )- Y( M) ) 

THD=TH (M )+RATD«( THN( K) “TH( M) ) 

EM3=‘XM3( .5 ,, 5, THD, THN( L) ) ' “ 

ydt=yd 

XD=( YN (L >~ Y( M) +SL0PE »X (M )-EM3*XN (L ))7( SLdPE-EM3) “ 

YD=YN( L) +EM3*( XD-XN( L) ) 

■"‘IF' (A8S‘(( YDT-YD )/'(YN( K) -Y Tm ) )“'). LT'. . Od'i)‘Gd' f O 23" 1 
KIP=KIP+1 

IF (K IP .G T. 20 > CALL E RR OR (2 30 ) 

CiO T 0 23 0 

23 1 RA TD =( YD -Y (M )) /( YN <K )- Y( M) ) 

2 23 5 cO NT iN UE 

OD =0 (M RA TD «( QN (K )- Q( M) 5 

TD=T (N )+RATD«( TN (K )”T( H) ) 

THD=TH(M)+RATDtt( THN(K>“TH(M) ) " 

PD=P(M)+RATD«(PN(K )-P(M) ) 

RD=R (M ) + RATD»( RN (K )-R{ N) ) ' — 

_ WM T= W( M) +Q (M )« Q{ M) /2 . 

T= WN'(K )+ (JN (K )« GN (K )/'2‘. 

WD T= WM T+ RA TD WK T- WH T) 

WD =W (M )+ RA TD «( WN (K )- W( M> ) ' 

IF(IeQ,Eq.i) Nd= wDT’“Qn«^D/2. 

■ ■ ■ DO 2 35 J =1 ,NSP 

__23 5 AL PD (J )= AL P( J, M) +R AT D*^ (A LPN( J, K) -A LP (J )) 

CA LL ■ A LL rcP XD , W D , RH OD V R D , GA MD ,E MD , X'MU DVQDVT D ;PD", A LP D )’ 

DO 1502 J=l. NSP 

ALPN (J ,L )= ALPD (J ) ■" — — 

150 2 ALPDUM (J ) = ALPN (J ,L ) 

IF(lEQ.eO.l) GO TO 3002 ‘ ” 

TERM2=RH0D<^0D 

i F ( B . G T : 0 7) TE RM 2= ( T ER M2 +R HO N ( L ) «0 N (T ) )/ 2, 

IF (B .ED, 0. )CPXN( L) =CPXD 

QN(L) = (PD-PN(L))/TERM2+QD ' “ 

TN (L ) = TD + ( ON (D + OD ) PN (L )-PO )«EI N/ (CPXD+CPXN (L) ) /TERM2 
CALL ALL ( CPXN (L ), WN (L )r RH ON (L )# RN (L ), GA MN (L ), EM'N( L) .X MUN( L) ,ON( L) 
_1_»TN( L> ,PN( L) , ALPDUM) 

" GO" To 30 03 ■■ ■■■ 

3002 CONT INUE 

RHON(L) = RHOD " «-'(PN(L)7PD ) « * ( 1 ,'7G AM D~ ■"') ^ 

_ VVl=QD ^OD _ +2-<^GaMD /(GAMD '1.)»(PD /RHOD -PN(1)/RH 

ION (L ) ) ' ■ " ' ■ ' 

QN (L ) = SQRT {_yVL > 

W N ( L' ) = W D " + . 5 » ("Q’D ' ’ VOD ' - V ’X ' )' ~ 

RN (L )=RD 

TN(L) = FT(PN(D ,RN(L) iWNCL) )■ • 


-AIA- 



GAMN(L > = FGaH(TN<U »PN(L) ,RN(L) ) 

" Bm ( L > =0N { L‘) /S'QR T ( GAMN ( ll )'t>PN ( L )‘/RHbN (T )1 
^ XMUN (L )= aT aN (1 ./ SORT (EHN (L . ) ) 

3 00 3“ CO"NT INUE “ 

IF{IdPT.EQ,0) 60 TO 8^ OO 

ine .GT70, )“G0 "TO 8401 “ " ' ~ ' 

A= ,5 

“e='.5 

GO T 0 25 3 \ _ 

8 40 0 continue 

IF (B .GT. 0. ) GO T 07 52 0 _ 

■“IFIL .NE. ISHOC) GO TO 6 06 

IF (L ,G T. 2) G 0 TO 6 50 _ __ _ 

"RAT= {YN( 1 )- YSHO C) /( YNT lb‘-Y'{M ))■ ~ 

y(m)=yshoc 

' xTmV =fsh oc ■' ' ' ' - ; 

OH 0= 0( H) 

Q "■ TM) =0' ■ “(M )+ Ra T* (Q “ N( T )-“Q‘ ■■ ■( M) ‘ " " 

T < M) =T (M )+ RA (T N( 1 )- I ( H) ) 

''R'" ( N)"=R‘"“ (M )+ RA T^»' {R ■ N( '1)-R ' ( M) ) ” '' 

CMl P aN DW (P (M ), pN (1 ), RA Tr P< M) ,W (M ), WN (i ), QH Q» QN (1 )» Q( M) ,W (M ), lE Q 

_ . . . . „ .. „ __ _ 

TH ( M) =TH (M )^-RAT« (Th N( 1 )-TH ( M) ) _ 

XMU( M) =XMU (M )+RAT^^ (XNUN( 1 V- XM U ( M) ^ ~ 

GAM(M) = 6AM{M)+RAT»(GAMN( 12”GAM(M)) ^ __ _ 

”■ DO'“655' J = 1,NSP - --- ---- - - 

_J>55 ALP{ J.M)=ALP< J,M )^-R AT<K A LPNiJ , 1 ) :lAIP > 

Go'foTob 

_650 continue _ „ 

RAT= (Y(f1-l)"YSHOC)/( Y(M-i)~Y (M) ) 

Y(M ) =YSH0C 

“ 'x(M)“=xsHoc' ■ ■ 

OHOfQ(M) _ _ 

— Q- • "( M) <M")+ RA TM (ir-rb- 0" “'( M)“) " ' 

T ^M> =T (M ')+RAT«(T (M -1 )- t ( M) )_ 

R • -( U)-rR- --- (|vj -)V^A T# (R‘" '(M "■ 

CAlE Pa^D>^ (P (N ), p( M- i_) ^PaT_,P (M )/r( m) ,W (^-l„)*_O^Hq^Q(_H- ), W( M) 

■ ■ " 1,1 EQTO ) ■■ 

TH ( M) =T H (M)+RAT<^(Th (M-I)-th ( M) ) 

XM7UTmT=XMU <My+”RA"tttl"XMU "(Tr-i ' ' ■“ - ‘ ~ 

gam(m)=gam(M )^rat« (Gam (M-i )^ga m( m) ) 

DO 6o5 "J=1,NSP ' " ■ “ ^ ■ ‘ 

605 ALP< J» M) =ALP (J ,M )+RAT« (Alp (j ,M-1 )- ALP( J# M) ) 

" 60 6 continue 

A= .5 

' "g - 5 ‘ - ~ - - - ~ 

GO TO 25 0 

7520' CONT INUE ' ” ' " " 7 

IF(IFENCE,EQ.O.OR,XN(1).NE.XcOWLH) go to 3108 

XXP= (YN( L) -BFENCE) /AFENCE+XBP 

IFCXNCL) ,LT,XXP) GO TO 520 

DUH=“rYN{X l-YNTky)/ ( XNrLT~XN'( K ) ) " “ “ 

XXP= ( AFENC£<^XBP-BFENCE-DUM«Xi\!(K)+YN(K_) ) / ( AFENCE-DUM ) 

3108 CONTINUE 7' ' ' • : - 

IF (X J1 ,EQ. 0. . and , I OVER .Eq. 0) GO TO 5 20 

IF (XN( L) .L T, XX J1 (NXX J> ,AND ,XN( L) .LT , XXP) G 0 TO “5 20 

IF (LST .GJ_, 0) G0__T0 5 20 _ 

lS'T=l'“ — 

XN N= XX P 

IF(XN(L) .G'E'. XX J1 (NXXJD" XNN= XX Jl (N'XX'J 





RAT= (XnN-Xn< K> ) / (XN( L) -Xn( K) ) 

IAA=IAA+1 

iFdpP. GT.f'iPTS. OR. IAA.lt. IPP)G 0 TO 3999 
WR It E( 6r 3'9 98 ) “ 

3998 FORHaT(» INDEXING IN_CHANG6 OF ORIGIN OV ER LA PS _1NI TI AL DATA lAA EO 

lUALS IPP<^) ■ ' ■“ 

STOP 

3999 XN(L" )=XNN " ' ' 

YN(L ) = YM (K )+RA T« (YN( L) -YN( K) ) 

TH’N(T r=THN (K ) + RAT« dHN'd )-THN( kT ) 

DUM1=AlOG( PN (K )) „ , 

DUM2=AL0G( PN (L ) ) 
dum=domi+r AT »( dum2«dumd 

■ PN(L ) = EXP(DUM) ■ - - 

IDPT=1 

'A=i. ™ ■ 

8=0. 

GO TO" 25 3 " ' “■ 

8A01 IDPT=0 

IF (I FENCE. EQ .0 .OR. XN (i )7NE .XCOWLH) GO' t0"9‘80 4 

IF (I aA ,LT. 3. OR .XF( 1) .EO. XFENCE ) GO TO 9804 

WR IT'E( 6/98 05 ) ' ‘ ’ ’ ^ ' 

9 80 5 format (« fence M aY N OT BE ENTIRELY SUPERSONIC - CHECK FLOW FIELD*) 

■ STOP ■ ■ ■ “ 

9 80 4 cO NT iN uE ^ _ 

X F( lA A) =X N (L ) “ ■ “ ' 

Y F( lAA) =Y N (L ) 

Q F( IAA) =Q"N (L ) ' " 

_R F( 1AA)=R N(L) 

W F( lAA) =W N (L ) ■ ‘ 

T F( IAA)=T N(L) 

P F(IAA)=P N(L) ■■ "" ”” 

THF( IA A) =T HN (L ) 

DO 8 40 2"’“J=1/NS P ■ 

8 40 2 AL PF (J , I aA )= aL PN (J ,L ) 

IF (I AA'.EQ, 2. AND. LS TT .E Q', 0)*G 0 'TtD 5 20 

LM AX =L 

60 TO '265 - 

_ 52_0 CO NT INUE^ _ 

■ I F' ( X N ( L )■ , L"r.'XX Ji'CN XX JlTA'NDTX N( n ':rT7XXPT 'G 0 ~T 0 ~ 5"52 0 

IF(LSTT^.EQ,0) GO TO 5520 _ 

GO TO 265 

6520 CO NT INUE 

A = l.-^ 

l=lmax 

KK=KMAX ■ ■ -- - 

K=LMAX-1 . . 

5 011 continue — 

IF (B .EQ. Q. ) THN( L)=T H( KK ) 

63 0 CONTINUE 

■ EMSL=XM3(,5. .5 ,TH( KK ) , THN( L)“) ■ 

EM2= XM2( A# THN( K) ,XMUN( K) ,ThN (L ), XMUN (L ) ) 

if(Idesgn;£o;o,and.iflip.n£,i)' - • 

ICALL BODL( X( KK ), Y( KK ), TH (KK) ,XN( K) ,YN( K) ,EM2 ,XN( L) ,YN( L) ,THN (L ) ) 

IF(lFLIP".NE,i) GO 'TO 631 ' 

XN (L ) = (YN( K) - Y (KK) +EMSL*X( KK )-EM2*XN (K ) ) / ( EMSL-EH2 ) ■ 

YN(L ) = Y(KK ) + EM5L*(XN(L )-X(KK )) 

-Al$~ 



__631_ CONTINUE __ _ _ _ „ 

Bl=Fl( AiB'* XMUNCK)7'GAMN(K) /PN(k)i XMUKUL) /GAHNa^^PNa) ) 

XDUMKf XN (K )-XOR( NX X^-1 ) 

XDUM‘L= XNin-XO'R(¥XXJ-l ) 

=F ?( A. -1 . , XJ ,X Jl ,X DU MK , Y N( K> ,T hN <K ) , XM UN (K ) , XD UMLr VN (L ) > TH N( U __ 

T/xhUni'd') . ■■■ " • * 

DU Ml =8 2-» (X N(_L> -X N( K) ) _ _ __ 

" 'IFCXUI .GT.“0. )DUMl=B2«AUOG( XUUML/XDUMK) 

^DUM=__ (THN (L )-THN( K) -D UMP /B 1 

PN TiI ) = PN (K )«EXP( dumT' 

_ „ IF (B .E Q. 0. _) GO J 0 18 

'Ll =1 MAX 

L2=KMAX _ „ _ — 

’ir (XN( Ll )TGT .X XJ 1 ("NX XJ )-T. E-OA ) GO T O 18 21 

Zl=AxX (NXX J-D ai X( L2 )-XBP) +B XX (NxX J~l) 

Z2 =A‘‘XX (N XX J- 1) «TXN (L 1) ~X'8P )+ BX X( NX XJ-^ 

DAX= (Z1+Z2 )/2. «ABS (YN( LI )-Y( L2 )) 

DAY= (Z1+Z2 )/2. *( XN (Ll) -X (L 2) ) 

PAV=<PN(L1)*^P<L2))/2. __ _ 

“PA ( P AV yp I N“-l . ) «P I N F« 14 4 , 

DX tH X= PA V» da X 

DYLFT="+PAV*DaY" 

XNS= (X N( LI X( L2 )) /2 ^ 1 - 

■“'YMS= (YN( LI )’^ Y('L2 )) /2 YSHF'T“ ' 

DM0M = YMS^»DXTHX-XMS^DYLFJ _ 

XTHX‘=XTHX + DXTHX ” 

YLFT=YLFT + DYLFT 

■“ XmOM = XMOM+M^^ 

1821 CONTINUE _ 

DO 1803 j=l,NSP 

ALPN(J,L) = ALP( J.KK) _ _ 

150y'ALPDUM ( J) = ALPN( J,L ) 

jFdeQ.EQ.D GO T0 3 00 4 ; 

' ‘ . TE RM 2= RH 0( KK > Q(' KK A 

1F(B.GT.0.) TERM2= TERM2+B«RH0N (L )i»QN (L ) 

ON (L ■)= (P (K K) »P N( L) )/ TE RM 2+ Q( KK") "" ^ 

IF(B.E0. 0. ) cP XN (L )= cP X( KK ) _ 

TN rC)=T( KK )“+ (QN( L) +G) (KK) ) »( PN (LT-P( KK )T5'E IN/( CP X“( KK ) + CPXN (L ) ) 7T 

ieRM2 , _ 

CALL‘ 'ArL“rCP'XNrL“)7NN (L )7 RHW{L )"/RN Vl ) 7GA MNTL )7TM>;i( L) /xHU NfL') ,Q'N('L)' 

1_,TN( L) ,PN( L) ^A]_PDUM]i __ 

” " GO TO 30 05 

3 00 4 CO NT IN UE _ _ 

RHON (L ) = RH0( KK )« (PN( L) /P (KK ))«<»( 1. /GAM{KK) ) 

VV L - 0 ( KK ) i^^( KK )^^2.»GAM(K K) /(GAN(KK)-l. ) P ( K K ) /R HQ ( K K ) -P N ( L ) / R H 

IonTl)) 

QN (L) = SORT (VVL ) 

■" WN (L ) = W( KK ) + ,5 «("Q( KK )oQ( Kk )- VVL) 

RN(L)=H(KK) _ ; 

IN (L > = FT (PN( L) »RN( L)7 WN( L )') ■■ '■ 

GAMN(L) = FGAM(TN(L) .PN(L) ,RN(L) ) 

FMN("l') =ON(L)~/SORT(GAMNTL )»PN(L )/RHONlLl )■" - • - - ■ ~ 

XMUN (L ) = ATAN (1 ./SORT (_EMN (L_)^»»2-l dj __ „ _ _ 

3005 CONTINUE ' " ■“ 

IF(IFLIP.NE.I) GO TO 632 _ _ „ _ _ 

■ P2=PN(L)/P(KK) 

CALL PM(P2,TH(KK )>EMC 1 , TH 2 T . EM C2 , ^ 1) _ 

ER4=TH2T-THN(L r 

IF(ABS(ER4).LT ,i,E-04) GO TO 632 

' IB = lB+i 



IF (! GT 45> Go TO 633 

■ IF (1 B. GT .2 ) GO TO 634 

£R6=ER4 

'THSVV=THN{L) 

THN< L) =THN (L )- ,0 1«AB S( THN( L) ) 

'go to 63 0 

633 WRITE(6» 635) . , .. 

63 5 P0RMAT<* error IN PM LOOP AT CONTACT IN MAlN«) 

STOP . 

63 4 DUHI =T HS VV “E R6 «( TH N( L ) ~T HS VV )/ (E R4 «E R6 ) 

ER6=ER4 

“THSVV=THN(L) 

THN(L)=DUM1 _ 

GO TO '63 0 

632 C0NTINUE__^ 

IF <B‘.GT^ 0. ) GO T01265 

“ B= .5 

GO TO 5011 

1265 IAAS = I AA ■■ 

IF ( I FL IP ,E Q. 2. AN D. XN (L ) . GE . X FI NAL)_J HA Lj_=l 

IF'VrHALT .'EQt 1) “G'0“TO‘ 160 2 

IF (XN( L) .GE. XF INAL ) GO T 0 18 02 

IF (X J1 .FQ. 0. .AND .1 OVER ,EO. 0) GO TO 2 65 

IF (1 fE nC E. EO .0 .0 R. XN (1 ). NE ,X cO WL H) GO TO 310 9 
IF (XN( L) .LT, XF ENCE > GO TO 26 5 
XN N= xF EN CE „ _ „ 

'GO TO 3110" 

3109 CONTINUE 

IF(XN(L) .LT.XXJKNXXJ) ) GO TO 265 

1802 CONTINUE _ 

XNN=XX J1 (NXX J) 

3110 continue _ __ _ 

RA T= (X NN -X (K K) )/ (X N( L) -X (O) ") 

IAA=1 __ _ — 

LSTT=1"'"- 
XF(1AA)=XNN 

* IF (1 FL IP .EQ. 1) GO TO 2 60 5 ' ‘ 

XCH=XNN -XINTL 

- LK=3 

IF (XCH .LT. XXL( 3) ) LK =2 _ 

IF (XCH .LT. XXL( 2) ) LK=1 

YF (I AA ) = AA1( LK ,2 )«XCH<»XCH^ AA2( LK ,2 >«XCH+ AA3( LK #2 ) 

THF( lA A) =A TAN( 2. «A A1 (L K, 2) «X CH +A A2 (L K, 2) ) 

' GO To 26 06 _ 

" '2 60 5 YF ( i AA )= Y( KK )+RA T« (YN( L)'-Y (KK) ) 

THF( lAA) =ATAN( TAN( TH (KK) ) + RAT« (T AN (THN (L ) ) -T AN (TH( KK )) 
2606 CONT INUE 

OF (I AA ) = Q( KK ) + RAT* (QN( L) -Q (KK) ) 

RF ( I AA ) = R( KK ) + RAT« (RN( L) -R (KK) ) 

PF ( I AA )= P( kK )+ RA (P N( L> -P (K k) ) 

Ll =l A A 

__ L2=KMAX 

Z1=AXX (NXX J-1) «( X( L2 )-XBP> +BXX (NXX J-1) 

Z2=AxX (NxX J-D «< xF (L 1) -XbP )+BXX( NX XJ -1 ) 

DAX= (Zl + Z2)/2. «ABS(YF(L1 >-Y(L2) ) 

DAY=(Zl+Z2)/2. o(XF(Ll)-Xa22)_ 

■ PAV=(PF{Lr)”^P(L2))/2. 

PAV= (P AV/P lN-1 . ) «P INF«144. 

DXTHX=PAV»DAX 
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f OYLFT=+PaV*DAY 

''■XMS=(XN( LI )-^X( L2 rr/2." -XShFT ' 

_YMS=_(YF( Ll )+ Y( L2*)) /2 YSHFJ 

DM OM = Y MS XT'HX -X MS eP YXTT” 

xthx=xthx+dxthx 

YLFT=YLF T+DYXr T 

xmom=xmom+dmom 

' WF (I AA ) = W( KK )+RATo (WN( L) -W (KK) ) " 

I TF(lAA) = T(KK)-^RAT<^aN(L). "T(KK) ) 

W A = W' < K K ) +0 '( K K ) « Q ( K K V/ 2 . ' 

IF (I EQ ,EQ. 1) I'iFC lA A) =WA-QF ( I AA { I AAy 2. 

'seAi j=i; NSP ■ " 

3e4l aL pF (j , I aA )= aL P( J» kK )+ RA (A lP N( J» L) -A LP {J fK K2J 

■■ IAA= I A AS ■ ■ ■■ ““ ■ 

__26 5 CONTINUE _ „ _____ 

'C‘ ■ ■■■■lF (P N( LM AXTTlXXPENXANDXTDE $S NTXQ ;i ) GO ‘'T O y~6 W" 

kn MA X _ _ 

BET=BETN 

EM Cl =E MC 2 

DO' 2 70 K=r,KHAX 

X< K) =X N( K) 

Q( K) =QN( K) 

■T(Kr=TN(K) ■' ■ ' 

TH(K)=THN(K) 

PVK‘)=PN(k) 

R_ _ N(K)_ 

Y “(K)=w“ ‘n"Ik) 

CPX(K) =CPXN(K) ^ 

RH0(K)=RHON(K) 

GAM(K) =GAMN(K) 

” ) =EM N( K ) ^ 

XMU^ Ki =XmUn( K) 

"W 2 71 J'-l »N.SP 

27 1 AL P< K> =A lP N( J» K)_ 

■ DO" 1 50 J= 1» NS P 

I 5 O 4 aLpDuM (J )=ALP( Ji K) _ _ _ _ __ •_ 

'CALL' all rCPX( K)",W (K )V RH 0rK)'7R (Kl^XA"Mrkr;E"M( K ,XMU <k )> (H K) ,T (K ), 

_1 P( K) , A LPJU M] _ . 

27 0' CONT INUE ” 

_7266 IF (N-NPTS) 7201, 72^^^^ 

GO To 50 0 

' ' 7202 CONT INUE “ 

_NSTAR=1 

”1PPSV=1PP+1 

_N = N+1 . _ 

■ GO TO 50 0 

9123 KTEST = KMAX/2 • ; 

KTEST1=( KMAX-1 )/2 

KT3 = KMAX-1 _ „ _ ' 

i F ( K TE S T :e Q ;k TES tD 'K T'3 ='KmTx" "■ ' 

1 = 0 _ ___ „ „ _ 

'DO 6A12 L=1»KT372' 

1=1 + 1 _ 

X (l)=x (L) 

Y__(I)fY (L) 

■ Q ( iT = Q (L )" "■ 

p(I)=P (L) _____ 

"T (I)=T (L> ■ 
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w (prW (L) _ _ 

p ( T ) =R (L > 

_th ( IJ =TH (L ) 

eM (ir=EM (L) 

CPX( I ) =CPX {L ) 

RHO( I) =RHO (L ) 

GAM( I) =G AH (L ) 

■ XMU( I)=XMU(L) 

__J)0 6414 J=1»NSP 

6414 'aLP< J* I)='ALP(J.L ) 

6412 CONTINUE 

IF (KTEST .EG.KTESTi ) GO TO 6413 

. I=I*1 

l=km ax 

X ( I) =X (L ) 

Y *n>'=V' (L ) 

Q (I)=0 (L) 

P ( I) =P ■ (L ) 

T ( I) =T (L ) 

W ( I) =W (L ) 

R n> =R (t- ) __ _ 

Th VI)=TH (L) 

( I) =E H (L ) 

■ CP X( I) =C PX (L ) 

RHO( D =RHO (L ) 

GAH( I) =G AM (L ) 

XMU( D =XMU (L ) 

' DO 6 41 5 J=lf NSP 

6415 ALP( J» I)=ALP(J,L) 

6413 KMAX=KMAX/2^1 
lmax=kmax 

60 TO 6060 

_fcN D 
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CO MM ON / C OW L/ I C OW L » MM , X M ( 9 ) ^ Y M ( 9 )_, P M ( 9 ) W M ( 9 ) » R M_M ( 9 J »_TH M ( 9 ) »_0 M_(_9 L» — 

_ ^ ^ ) 7gM ('9 ) » XM'UM ( 9 ^ GH M ( 9 > , A LP H ( 7 » 9 ) » C PX M '( 9 ) . , , 

COMM ON/SHAPE/ Al(3»2)» A2(3_f2 ),A3(3 r2)» ^ X_U^_3 XL X 

COMMON/L TM/XSHFT',YSHrT ,XTH'X» yCFT , X MO M 

CO MM ON /X XJ /N XX J1 , X XJ 1 ( 6_) f A X^( 6 ) / BX X ( 6 ) i X OF^ ( 6 ) 

^0 mM OM /X Fl NA L/ Xf* IN AL ' ' \ ’ 

COMMON_/_SFV NSP _ 

cOm^oN/^/ fi NTcPiNtFio 

CO MM ON /B/ WTMOLE 

qOMMoN/D/ GA Ml NF ,E Ml NF ,R INF » wI NP 

C0MM0N/F/A9, 69,1 800, XWFfNBODfYEND . - • . Kiw vdp 

C0MM0N/ETX/XJ,XJ1, NPTS ,10, IR EFL, IC hE M, IP UN CH , 1 DE SG N. ^ r ^ , 

1 Y6 P, TH BP ,R AD ,X BO D, YB OD ,T HB OD ,R AD B, XE ND ,T HE ND »R IH ,_Y E.X ll .T hSJ_, TE_S^_ 

llRFnY07RADB2;RR AD (2 0r; NS TA R.'YN OZ ,E IN , ^ I' , 

CO MM ON /X F/ XF (2 00 ) » YF (2 00 ) 00 ) » OF (2 0^_) » IF_1.?._^.P li- 'Ij fJ 

^ CO MM ON /X / X( 20 0) , Y (2 00 ) » P( 2^ 0)_»Q (2J30J_,_T^ 

CO MM ON /X cO Is'L /X CO WL 

MM ON /F V ^ ^ O0')7EMrr20'Q')TX MU'F 720 O'FTFF f215'0 r, RT'r20 0 )TO AH F ( 20 0 ) » 

2XMASSF(200), H F ^ A ^ A?- AL-lAi-L.PlA?.2 ' 

3THETA r20 )■"■ 

COMMON/VAR/RHO (200 ) » __ , . - . 

VpM(200 r,'XMU(200 )^ CPX(^2 00 ) , H ( 20 0 ) , R ( 200 ) . GAM ( 20 0 ) ,XMA5S(200) , 

2 XN(200)»YN(200),QN(20 0 ),TN(2 0_0)_»PN(2_0 0),ThNj(20 0 ) R HO N _ 

3T2”o6T. 1TTnT 260 r. XHU¥( 20 0 rrCPXN C 20 0 )7^N"( 20 o') ,”R¥ (200) ' , on , 

4 X M A S S N ( 2 0 0 ) , A L P N < 1 0 . 2 0 0 ) , S H 1 0 ) . . H 1 ( 1PJ_! E 1 ?-? ’-'i ' 

5,ALPDUM(10) 

C0MM0N/1E0/1EQ»P1N»R H OJ N F/_y ] N F »_P I N F 

‘ “common "/IOVER/ IOVER 

CO MM ON /X TJ 1/ XT Jl — 

CO mM ON Vr EN CE’/l FE NC E^ AF'EN CE ,B FE NC E, XF EN CE \ -r^/ nr- , .1 n 

DI MENS ION WTMOLE { 1 0 ) , T YP E( 6) , T YP El ( 4 ) , TYPEJJJ^) ,J YP E;5 ^ 

OA TA T yPE2 '7iOH ' N07ZLE‘~ f IOhCENTERBOD Y/ 

data T YP El /I oh HY dr OG EN A , 1 QH HY DR QC AR BO , 2 HI R , 5H N AIR/ 

- OATA 'TYPE/IOHTWO D I M EN S , 10 HA X I SY MM ET R , 10 hLTN E 50 UR C , 5H IONA L , 5H IC 

'“O^TA^TTPE57rOHFnL TlOH D'E'~S »9HF I E L 0,9Hf G N 

data T YP E4 /I oh FR 07 EN , 1 OH EG U I L I BR I U » 1 RM.X 

data XXU'/3«1 .E+0 6/ ,XXL73ttl ,E+06/ »EN0 7l .E + 0 6/ 

data IOVER/0/ ^ ^ 

400 F0RMAT(1H1 ) ■ 

16 F 0 RM AT ( 3 X , M , 3 X , 1 0_E1 2 ,_4 ) 

“T602”F0R^A^'(i0X ,loE"l2 ,4') ‘ 

lF<10VER.£Qr2) 00 JO 916 

10=0 

irefl=o _ _ 

IPUNCH=0 

_IDESGN_=(^ 

iCHEM=0‘ 

RE AD ( 5 , 6 89 5 )_ , J 1 . J2 » N PT S / I E Q , I C OW L i- 1 0 VE R 

5100 NSP=7 

6 89 5 FORMAT (161 5) 

■■■ ' RE AD (5 ,6 3) XBP ,X80U, XCOWL, RTR, 1BST ,Xf I NAL, XT Jl 

RE AD ( 5 , 6 3 ) X SH FT , Y SH F T > X TH X , Y L FT ,XMQ M 

'READcS^^a^^) NXXJI"' 

DO 9393 I = 1»NXXJ1 _ . . 

— 9393 READ (5 f 9463) XXJK I) 1 AXX ( I ) , 3XX( 1) ,X0R ( I ) 
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9463 format (4E10. 0) 

READ (5,4)' IF EN C£ , A FE NC BF ENCE .X FE NC'E 

A F0RMAT(I5.5X3F. 10.0) _ 

XX Jl (N XX ji ♦! )= XF IN aL ' " 

AXX( NXXJl+1) =AXX (NXX Jl) 

BXX( NXXJl+1) =BXX (NXX Jl ) - ■ 

XOR (NXX Jl+1 )=XOR (NXXJl) 

"' ■■iFCIFENCE.EO.O) GO TO 321 

XX Jl (_N XX Jl 4*2 ) = XF INAL 

AXX (NXX JIV 2 ) = AXX (NXXJT)' " ' 

BXX (NXX Jl+2 ) = BXX (NXXJl) 

XOR (NXXJ1+2) = X0R (NXXJl) ' " 

321 CONTINUE 

XI NTU=XBP '■ 

XINTL=XBCD _ 

READ (5 /6895') NUWSE gTNL WSEG 

READ(5f6 36 3) (XXU(L) fAl(L#i) ,A2(L»1) »A3(L,1) »L=1»NUWSEG) 

RE AD (5, 6 56 3) ( XX L ( L ) , A 1 ( U , 2 ) , A 2 (L , 2 ) . A 3 ( L , 2 ) , L =1 r N LW SE G ) 
6 363 format (4E10. 0) 


3531 WRITE( 6, 7329 ) TYPE5( IDESGN+1 ) , TYPE5( ,IDE5GN+3 )' ‘ 

2329 J0RMaT( ^ 41X*N_ 0 J_Z__LJE -__CENTERBODY « // 50 XA IQ » A 9) 

JPl= ji+1 +2<^J2 ■ ‘ “ ' 

JP 2= JP 1^ 3 


W'R IT E( 6f 7330 ) 

1 ,TYPE1 (I CHEM+3 ) , NP TS ,RTH 


TYPE(JP1 )»TYPE”(JP2 )'»TyPE 1 (ICHEM+l) 



7 33 0 FORMAT! /3 7X.»F 0 R“ N 0 N U'N I F O R M GAS' F L '0 W 

2»/// ^ 10X.«TYPE OF FLOW IS ^»,A10,A5. ^ FOR 

3 A1 Of A'5/1 OX ,«NUMB6R OF POINTS 0 N IN I T 1 A'L DA TA L INE ' IS , I 3/j 0 X,' ' 

4«THR0AT RADIUS (RTH) = -».E13.5) 

WR ITE( 6» 5001 ) TYPE4( lEQ^-1 ) , TYPE4 ( IEQ-^3 ) ' ' 

5001 FORMAT! 


IIOX^CHEM ISTRY is «A1'0,A1) " - - . ~ 

WR 1TE( 6f 7500 ) XC0WL,XF1NAL 

7 5 O 0 FO RM aT (i oX »C oW L 'jR Al Ll NG E'dG \S « El~3'/57i0 X* AXTA‘l CO OR DI NATE 0 F' EN 

ID OF RUN IS «E 13 .5 ) 

IF (J 2 , EQ .0 ) GO To 7510 "T - 

WR ITE( 6, 7501 ) XTJl 

750 1 FORMAT (1 OX «A XI AL COORDINATE OF S TA RT 0 F CA RT ES I A N $Y ST £M IS «El3,5 
1 ) 

WR IT E! 6/75 02 ) ' 

7502 format (//25X^»L INE SOURCE COORDINATES*) 

WR ITE( 6f ^505 ) - - 

DO 7 50 3 1= 1, NXXJl 

7503 WR ITE( 6, 75 04 ) XX Jl ( I ) , XX Jl ( Ui ) , AX X( I ) ,B XX ( I ) , XO R{ I ) ' — 

,J504 F0RMAT(10X,4E11.3f 5X,E11.3) 

7505 FORMAT ( 1 5X <^X «4 X* TO *4 X* X* 11 X«CO ORD I NA TES*'l‘?X* OR 1 G I N' XV) 

7510 CONTINUE 

WR ITE( 6, 2020 ) " ‘ ‘ ■ ' - '■ 

2020 FORMAT (//25X«UPPER WALL COORDINATES*) 

WR ITE( 6* 20 71 ) - - - - • - 

2 07 1 FO RM AT ( 1 5X *X »4 X» TO *4 X* X* 15 X* CO OR Dl NA T£ S* ) 

20*^2 form aT (1 OX ,5e1 1 / 3) : 

DO 2 05 0 I = liNUW5F.G 

IF (I ,EO, 3) Go TO 2 051 ‘ - - _ 

WR ITE( 6 , 20 42) XX U( I ) ,X XU ( H-1 ) , Al ( I ,i ) , A2 ( I , 1 ) , A3 ( I ,1 ) 

GO To 20 50 ‘ - - - 

2051 WRITE(6, 2042 ) XXU( I ) ,END , Al ( I , 1 ) . A2 (I , 1 ) , A3 ( I , 1 ) 

2050 CONT INU6 - 

WRITEC6, 2010) 

2010 FORHaT (//25X.L0WSR WaLL CO OR Dl NA TF S»‘) " ' ‘ ■“ 
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_ _ WR ITE( 6» 2071 ) __ 

■ bO‘'2 04o"' l = l* NrwSEG ■ 

IF ( I .EQj^3) GO__TO 2041 

‘"WR ITe( 6^ 20"42 ) XX Lrr> ,X XL '(T-^f) ,A1 ( 1 , 2 ) i T2 ( I , 2 ) , A3 ( I , 2 > 

_GO TO 20 40 

20^1 WR ITe( 6. 20'42~) yriTD ,END , Al (i ,2 ■)» A2 (I ,2 )i a3 (I ,2 ) 

204 0 CONTINUE ...... 

dO 7373 l=i» nuwseg 

J37 3 XXU( L) =XXU a )-Xl NTU 

' b’O’ 7 37 4‘I=1 >'nLwSHG ' 

7 37 4 XXL( L) =XXL (L )-Xl NTL 

3 53 2 CO NT INUE’ 

RE AD ( 5 , 6 3 ) E M I NF , Tj N_,_W j N F , P i NF 

DO 5 oO 2 i= 1 * NP TS 

P0_5 00 2 NS P 

■‘5 001 'aL PFl 

READ (5 ,6 364) ( XF__( I_ ) , YF ( I >,_P_f;n ) HJ 

IF (I'eQ .EQ. 0) 

IRE AD (5 ,6 36 4) ( (ALP F( J, I) ,J_=1 ,NSP >, 1,_NPJS 

6364 format (7E10, 0) ■■ 

IF (I COWL .EO. 1) READ (5,8698) MM, PM (MM) ,QM( MM ), TM (MM) , 

rt H M ( HM )' , RM (Mm) 

IF (I EQ ,NE. 0. OR ,I COWL .NE._1J _GO TO 7 511 

' WR IT E( 6/7512 ) " 

7512 F0RMAT(<^ FROZEN EXTERNAL CASE NO T ALLO WED* /■» CHE CK_ I NPUT MANUA L_FO 

'IR POSSIBLE USES OF Te'Q 'AND TcoWL^t) ' 

_ STOP ^ 

'l5Ti~~C0NT l NUE’ 

8696 FORMAK I5,5X,5E10 . 0) _ _ 

■ IF ( I COWLVEQ, 1) WR ITE(6750 05 )MM, PM (MM) , QM( MM ) , TM (MM) ,THM (MM) , RM( MM ) 

5005 FORMAT (//20X*DATA AT C OW L« /I OX *N UM BE R OF POINTS IN PRA NDTL -MEY ER F 

IAN is *I 2 /inx*PRESSURE‘ =^^E13: 5 / 10 X* VELO Ct TY =*E 

113 ,5/1 OX t^T EM PE RA TU RE _ _ = ^ ^ i?_^ ™ ^ 

~i--^O X* FU EL'“'Tb“.All RA Tro"=* El37^^^^^ 

IF (I RE FL ,E 0. 0) GO T0__4 00 0__ , __ 

““916 CONT INUE 

DO 4 00 5 I-1,NPTS_ 

lT=NpTS'l+l“ " " *' ' ' . 

x( II - - 

YnT) = -YF(T) 

_ 0( II ) = QP ( I ) 

f( II )=TF ( ! )“■ 

_P(n) = PF(I) 

TH ri 1) =-THF( I ) 

__R7in = Rr (j.) 

W ( 1 1 ) = WF ( i ) 

D0 4010‘J=1,NSP 

'4010 ALP( J, II ) = ALPF (J,I ) 

4 00 5 CONTINUE ' _ 

DO 4 00 6 I='l,NPTS 

XF (I )=X( I) __ 

YFd) = Y(i) 

0F(I) = Q(I) _ _ 

TF(I ) = T( I)"'”^ 

PF(l ) = P( I) _ _ _ _ _ 

■■■ THF( I)=TH( 1) 

RF(1)_=R(I) ... _ 

■ i') = w(i'>" 

DO 4 00 7 J=l. NSP _ • _ 

'4007 ALPP(J» I )=ALP( J, I ) 
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4 006 CONTINUE _ _ . .. 

DUM= -YrP 

_ . ybp= -y bod 

“ ' Y8 OD =DUM 

DUM=-THBP ___ _ _ _ 

'■ thbp=-thbod 

THBOD=DUM _ _ _ 

' 'Vend =- yend ' ' ' ' 

TH EN D = - T HE ND ^ 

iF (^ QVeR .EQ» 2) ReTuRn 

4 00 0 CO NT INUE _ 

XJ=Ol 

_ _ XJi=J2 __ „ 

■■ WTMOlE (1 ) = 1. 00 8 " ' ‘ ' ’ ' -- - -- — 

WT MOLE (2 )=16 , 

"wTMOLE (3 ) = is;016‘ "■ 

WTMOLE (4 ) = 2. 016 

WT MO lE (5 32 . 0 

WTMOLE (6 )= 17 ,0 08 _ _ 

WTMOlE (7 )=28 .014 ‘ ' ' ' ' ' 

WTMOLE (8 )=44 .011 

"wTMOlE (9')='28 .011 ' ■ “ ' ‘ 

WTMOLE (10) =4 4. 1 

R0 =1 .9 87 * ■ ■ ■ 

CA LL C OE FF (5 ,T IN , A Z, B2 iC Z, DZ ,£ Z, FZ ,G Z) 

CP =< AZ +8 Tl N+ cz *T \H p* \)Z *T iN 3"^ EZ »T jN 4) 0/ WT MO LE (5 ) ‘ 

_ CALL COEFF (7 ,T IN ,AZ» BZ ,CZ, DZ ,EZ, FZ ,62) 

' CP 1 1 = ( AZ +8 2«TI N^'CZ^T IN»*2+PZttT IN «« 3+ EZ -»T IN 4 ) 0/ WT MO LE'(7 ) 

RINF=R0/W1NF 

CPlN = . 232*CPIN-*-.768»CPII 

GAMlNFrl ./(i .-RINF/CPIN) 

EIN= (GAMINF-1.) EMINF<m» 2 
PI N= 1. /GAM INF/ EM INFO »2 
WRIT E( 6» 68 99 ) £M IN F, Tl N, 'Wl NF ,P IN F’“ 

6 89 9 FORMAT (// 5 OX IN FI NI TY C ON DI TI ON S» /5 OX ^ / 

1/40X,oMACH ,E 13 .5/4 OX TEMPERATURE (DEGREES 

2 K) ,E13 .5/40X ,<jMOLECULAR WEIGHT ,E 13 .5 /4 OX 

1<»PRESSURE (PSI) «E13.5) 

WRITE(6,6723) 

6 723’ FORM at (///40X, ^OUTPUT VARIABLES ARE* /4 OX , NO NDl M EN SI ON ALI Z ED ^^/ 

140X,ttAS FOLLOWS -»//40X,«X By RTH»/40X,«Y BY RTH«/40X,»Q BY FREE S 
2TREAM VELOCI TY«/40X, i^T BY FREE STREAM TE MP ER AT UR E« /4 OX , « P BY FREE 
3STREAM PRESSURE^^) 

IF(IEQ.EQ.O) WRITE{6,5003) ' " ' ' ' 

_5003 FORMAK /40X,*ALP(J) IS MASS FRACTION OF SPECIES ) 

DO 6897 1 = i7nPtS'"“ '' 

PF ( I ) = PF ( I )<^PI N _ _ ' _ 

'6897 CONTINUE ■" 

UI NF =1716. IN^l ,8/P IN 

RH0INF=(PINF /TIN/1, 6/1716, )»2116. 

IF (I EQ .EQ. 0) GO TO 293U 

DO 2929 I = lfNPTS ■ ' ' ' 

2929 WF ( I ) = FH (PF( I ) ,RF( n ,TF( I ) ) 

"2930 continue ■' ' ■* ' 

Hi 6= GAM I NF/W INF/ 2. *FM1 NF i»EMI NF 

H17=1, /H16 ' ' ■ 

^63 FORMAT(8E10, 0) 

IV'FORMaT (7E1 1. 4)’ '■ 

RETURN 

END ■■ ■ : ■ ■■ ■ ■" 
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SUBROUTINE C 06 FF ( I J , A ,_B , C D _ 6 , F G _ \ 

’''TMslOOO ” ‘ ' 

20 

10 GO TO" '(l5ri6,13Mlii2, 17 ,14i i'8 ,19) ,I 

_1_1 A = 2. 84 60 84 9E 00 _ 

'B “ 2 '4. 19 3211 6E “0 3 

C _ =- 9 . 61 19 33 2E “0 6 _ 

D ="9‘, 5122 66 2E-0 9 

=-3._^93 42iE_-12 

F =-‘9. 672537 2E 02 

G_ r-1. 4117 85 0E_0 0 , 

GO to 40 

12 A _ = 3. 71 89 94 6E 0 0 

B’” ‘V-2. 5167 28 8E-0 3 ' 

C =_6. 5837353E-06 

D =-8r29 98 7l6E“09~ 

6_ = 2. 70 8218 OE “12 „ 

p =-i: 05 76 70 6E 0 3 

G = 3. 90 80 70 4E 00 _ _ 

GO TO 40 ■ 

1 3. A 4 , 15 65 01 6E 0 0 • 

C =_ 5. 69 82 31 6E-0 6 ' _ _ 

D =- 4. 59 30'04 4E“0 9 

E 1. 42 33 65 4E^12 

” ' p =- 3^ 02 88 77 OE 0 4' 

„G ^ 86i6_^6E^0_1_____ 

“GO’ TO’ '4 0 " 

14 A _= 3, 6916148E 00 

=“1'.'3332552E-03 

__ _ C = 2.6503100E“D6 _ _ 

’ ■ 'D =- 9, 76 88 34 IE -10 

e =-9» 9772234E-14 ; 

■— -y~ =-T7’06 28 33 6E 0 3 

G _ = 2. 26 74 98 OE 0 0 

GO TO 40 “ " 

1^5 A = 2. 50 OO OO oE 0 0 

“”B =’“o7o’" 

C = 0. 0 . 

■ - D = ’ O '. 0 

6 =_0.0 

F' ”= '2T5'4 70'49 7E 0 4”' “ 

G_ =“4, 6001096E“01 _ 

60''TO '40 ’ 

1_6_A - 3, 02 18 89 4E 0 0 

8“ =- 2. 17 37 24 96-0 3* 

C = 3. 7542203E-06 __ 

0 E- 2/99 47 20 06 -0 9 

E = 9. 077754 7E-13 __ • _ _ 

r = 2.913719 06 0 4 ’ 


G _ = 2. 64 60 07 6E 0 0 
“GO t o “40 


17 A 

= 3. 82 34 70 8E 0 0 

B 

=-l. 11 87 22 9E -0 3 

C 

= 1. 246681 9E-06 

'D 

• = -27.10 35 89 66-10 

E 

=-5,25465516-14 

_ - 4 ' 1 ^ cr o ^ It T I- “ r» “7 


^ 3:58 52 78 76 0 3 


G = 5. 8253029E-01 
GO TO 40 ■ ■ 
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16 A=2* l7 0l 

P = i:0378U5E-0? 

C=-l .0733938E“05 

D=6 . 34 59175E-0 9 

E=-l ,6280701E-12 _ _ 

““ F= «4 .8352602E+04' 

G= 1. 0^6"^ 388 E +0l 
’ GO TO 40 

19 A=3. 78 7133 2 ' 

B=-2 .1 70 95 26E- 03 

C-5. 0757337 E-O 6 

D=-3 .4 73 77 26E-09 

E=7. 72i6 84iE-13 

F=-l .4 36 35 08E+04 

G= 2. 63 35 45 ^ 

GO TO 40 

20 GO TO (25, 26 ,23, 21 »22f 27 ,24, 28 ,29) , I 

21 A ■ ' = 3 . 04 36 89 7E 0 0 
B = 6. 11 8711 OE -0 4 

C " =- 7 . 39 93 55 IE -0 9 

D =-2, 033190 7E-11 
E" =”2 745 95 79 lE -1 5 

F =-8.54 9100 2E 0 2 

' G ■ =-l. 64 8133 9E 0 0 
GO TO 40 

22 A = 3.5976129E 00 

B = 7.8145603E-04 

■ ~ "C =- 2;2386670E-07' 

D = 4.2490159E-11 

E =~3. 3460204E-15 ” 

F =-l.l927918E 03 

G = 3.7492659E 00 

GO_ T 0 40 

2 3 'A ~ =^ 2.'67 07''53 2E‘T0 

8 =3. 0317U5E-.0 3 

“C ' '=- 8. 53 5157 0E-0 7 ' • 

D _ = 1- 17 90 85 3E-10 

E • ‘V-6. 19 73 56 8E -15 

F =-2- 9888994 E 0 4 

G = 676838391E' 00 ■ 

GO TO 40 

24”A " = 2.3545761E 00 ' 

B = 1, 5976316E-03 

C ’"=-6. 2566254E-07 

™P_ _=_ 1. 1315 84 9E-10 

E =-7. 68 97 07 OE -15 

F =-8;9017445E-4-02 

G = 6.3902879E 00 ' ' ' 

GO To 40 

25 A = 2.5000000E 00 ' ' ' 

B = 0.0 

c = oVo" 

D = 0. 0 _ 

E =0.0 

F = 2.5470497E 04 

G - = -4.6001096E-0l ' 

GO T0_40 

2 6' A " =■ ■ 2.53 72 56 7E'T0 

C _=-8. 80 17 92 lE -0 9 

D = 5. 96 43 62 IE -12 ~ 
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_ E =-5-57 43 60 qB 6 
F "■= 2.‘92 30 00 7E OA 

G = _4 . ^167_94 2&_ 0 0 

GO To' 40 

27 a _ = 2. 6895544E 00 
■ B ■" '= 9 , 98 35 06 IE -0 4 

C _=-2. 1879904E-07 

p- =“ i'. 98 02 78 5E -l‘l’ 

84 52 94 OE -1 6 

F = "3 . 88 li 79 2E 0 3 
___G = 5. 5597Q16E JO 

G0"fo“'4 0 

28 A= 4, 41 29266 

B= 3.19 22 89 6E -0 3 ' ‘ 

C=*l_ .297823E-06 

J=“2“. 41 47 44 6E ■-'10 
_ £= ,6 74 29 86E- l4 

F= “4 .8 94 40 43E+04 
_G= -7 .2 67 57 6?E- Ql „ 
GO TO 40 

9511519 

0= iT55 25 57‘67 E~ 03 
.19114 11E-Q7 
J=l. 13 50 33 6E-10 
E=-7 .788 2732E-15 

“ "'F= “r.4 2318 27E-04 

G= 6 . 53 14J5_ 

'4 f RETURN 

END , 



'SUBROUTINE ERROR(I IIT) 

WR IT E( 6, 10 UJ i i II 

io o' format (7 hi ER RO R= 1 5 ) 

__ CALL EXl T 

■return “ 

END 




'sUbROUTI NE’ BoDl('xI",V I* TH I, XL ,YL. EM2, xn 

COM M p N / S H A P E / A i_( 3 »_2 ) i.A 2 ) » A5 (_3 » 2 ) » 

Xr=Xt-xiNTL 

XL=XL-X1NTL 

' ICK=1 

ASL=TAN(THI ) 

XN = ( YI_-Y L+ XL t»E M2 -X AS L ) /< EM AS_U 

IPYxn'; lt ,XXL (3 ) ) L=2 

IF(XN,LT.XXL(2)) L=i _ _ 

YN =A 1< L » 2) <^X N« XN + A 2( L# 2) *X N-«- a3 ^ L » 2 ) 

_ YT =YL+ (XN-XL )»EM2 _ _ 

■'er = ( yn-yt) 

IF (A BS (E RKL T. GO J O Ip 

lCK=ICi<+i 

_ IF ( 1 CK ,LT, 15 ) GO TO 3511 _ _ 

WR ITE( 3512 ) 

3 512 rORMAT(« TOO MANY ITERATIONS IN BODL«') 

STOP ■■ 

1511 J F ( I CK .GT- 2 ) G 0 TO 3 513 

ER1= ER' ' ' '' 

XSL1=XN _ 

XN ’ '=1*"01*Xn ' 

L=3 _ _ j 

“■“go T 0“ 5 “ 

3 51 3 DU M= XS LI -£ R1 «( X N - XS U. )_/_( E R- ERJJ 

ER1= ER ” 

XSL1 = XN 

“XN =DUM 

L = 3 _ _ - 

GO TO 5 

10 TH N = A T A N ( 2 , » L_,_2 )<>XN '^A 2 ( L> 2) } 

■ “'xr=xY+xi nt‘l ' - 

XL=XL+XINTL ^ 

XN =X N+ Xl NT L ' " 

RETURN 



SURRoUTIwE BoDlK Xi yTtH )■ " ' 

COMMON/SHAPE/_ A1 (3 #2 ) r ^ * 2 ) » XXU( 3) #XXL )i XI NTU, XI MTL 

X=X-XINTU ■ - - - - . - - 

L-3 _ 

IF (X ,L T. XXU( 3) ) L= 2 

IF (X .LT, XXU( 2) ) L=1 _ _ 

Y= Al (L ,1 )^X*X+ A2 (L ,1 )«X'^ A3 (L ,1 ) “ 

TH = A T A N ( 2 . «A 1 ( L » 1 ) »X +A 2 { L r 1 ) ) 

x=x+'xintu' ’ 

RETURN 

E^^D■ “ ' ‘ ‘ ' 





SU qR 0^ T I nE C oW l < QE t ) 

CO MM ON /C OW_L/_I C 0 W L # _MM , K M (_9 ) » Y M ( 9 ) , P M ( .? M 1.?. > .* R HM ( 9 

1 RM ( 9 ) f TM (9 ) » GM ’( 9 ) ,' XM" UM (9 )7eM M (’ 9 )7 A lP M <"? T9 ) ’ i C PX M ( 9 ) 

C0MM0N/lEQ/ie0,PIN,RH01Nr,UlNF»PlNF_ __ 

"'common /VaR/RhO (2 00 ) I ' ' ‘ 

lEM (2 00 ), XMU( 20 0) ,CPX (2 00 ), W( 20 0) ,R (2 00 )» GAM( 20 0) ,XMASS (2 00 ). „ . 

■ "a XN (2 00 >/ YN (2 00 ) I QN (2 OO ). TN (2 QO )» PN (2 QO )» TH N( 20 0) ,RH0N 

3J2 00 ), EM_N( 20'0) ,XMUN(20O ,CPXN( 20 0) , W N(_2QJ ^RJ1(_20J_)_.G AM N( 20 0_) . 

" 4xM’ASSN”(2oO)7ALPN‘(10"/20 0r,Sirio7M^r(10’> ~ PR R( 20 ) . TE MP (2 0) 

__ _ _5 ,ALPDUM( 10 ) _ 

cOMM on'7 sp7”nsp 

COMMON /X/ X( 20 0) lY <2 00 ) P£2_0 0) ,0 (2 00 0)_ J_H( 20^^^^^ 

" XM (M M) =X (1 )" ' 

YM (M M) =Y (1 ) _ 

7m (M M) =P M( MM7-» PI N 

WM (MM) =F H( PM (MM) ,RM(_MM ), TM CMM) ) 

~D0 4 38 5 J= 1» NS P 

4385 ALPDUM (J )=0. _ . 

“call all (QQ» WM (MH)'iR'hM (M M) »RMi: MM ),'gM (MM) ,EMM (MM) ,X MUM( MM ), QM (MM) , 

ITM (MM) ,PM( MM )» aLPOUM ) __ 

BE T=‘(T HMYMM) -O'PT oXMUMC MM1+ TH 

ITT=1 

IFaN=M'M-3“ 

3JT=1 _ 

ii=MM-r 

VI =QM( MM )<^COS( BET- THM( MM )) 

U1 =QM(>.M')^ SI N ( BE T-THM(’MMT) 

„ _ U1 = ABS(U1) _ 

’GM1=GM (MM) -1 . 

GP1 = GM(MM)^1. 

XM1=EMM( MM )#Sl N( BET-THM( MM ) ) 

XMS=RHM( MMI )J»U1 _ 

IF ( rt ♦ EQ’ .1') ' UM =ll IM G Ml Ml h 1 +2 , ) /G PI /X Ml Ml 

7 RH2=XMS/UM ; 

''>M (I 1 ) =X MS ■» ( U1 -U M ) +P M ( MM ) 

VI =Q M( MM )» QM (M M) ‘ ' 

V2 =V T« VT +U M^ UM 

WM ( 1 I ) =W M ( MMO+JV 1- V2 )/ ^ 

~ ‘ RM'd i ) VR M'( MM )■ 

_RHM( II ) = RHEO(WM( II ) , PM (_I I) ,RM(JI ) rJM ( I I)J^ 

=RH2-RHM (I I) 

__ IF (ASS CER) ,LT. 1. E-04 )_^GO_TO 9 __ ^ 

'Tt = it+i 

IF (I T, GT .10) GO TO 100 

lF‘(TtVGt;2rGO TO' 'll 

ER2 = ER _ 

■“U2=UM 

UM=UM«.99 _ _ 

■ GO TO 7 " 

100 WRITE( 6,20 0) _ _ 

— 20 0^formaT (^ error ' in' hUg'on1(dt Loop in cowl») 

STOP . 

il DUM2 =U2“ER2^ (UM-U2')/ (ER-ER2) 

ER2=£R _ 

'j2i=U'M™"" 

_UM=DUM2_ 

go’t'o"?"" 

9 THM( II ) = BET-ATAN (UM/VT )_ _ 

'■ qH (1 1) =S QR TCV2) '■ ' ■ ' 
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XM ( I I ) =XM( MM ) 

YM ( I I ) =YM(MM> ■ ■ 

GM n i ) =h GAM( TM ( II ) » PM ( i 1 ) , KM ( I I) > • 

EMM(’il )'=OM(n)VSQRTrGM(ri'nrPM( II‘)/RHM( li )T ‘ 

XMUM ( I I ) =ATAN( 1. /SQR T( EHM( II )« «2 -1 , > ) 

THM( II -1 )=THM( II ) -- - -- 

CM (I I-l) =QM( II ) 

pM (I I- 1) =PM( ll ) ‘ 

WM (II-l) =WM< II ) ^ 

RHM'( Il -l ) = RHM( II ) 

RM(I I-1)=RM( II ) 

'tm (I i-i) =tm( 1 1 r “ ' ' 

XM (I I-l) =XM( II ) 

yM ( I 1“ !> =Y M( 1 1 ) 

GM(I I-l) =GM( II ) 

' pM M( il -1 )= ll ) ' 

XMUM (I I- 1) =XMUM( II ) 

OP =( p( 1) -P M( II )) /F lO aT (I fa N- 1) 

11=1 _ 

■ X M (I I) =x a r‘ 

Y M (I I) =Y (1 ) 

■ w M (I ir=w" (1 r 
p M ( I I) =P (1 ) 

■■ RH M(1I)=RH0(1) ' " ~ 

G M(1I)=GAM(1) 

Q M{ I I) =Q (1 ) 

TH M(II)=TH (1) _ 

XMUM( I I )=XMU(1) 

T M(II)=T (1) 

R M(I1)=R (1) 

EM M(II)=EM (1) 

DO 12 LL=2fIFAN 

N=LL 

"* kK -1 

XM (N )= XM (K K) 

ym (N )= ym (Kk> " 

PM (M ) = PM (KK) -DP 

' AL NR =A lOG( PM (N )/pM (K K) ) ' ’ “ “ 

RH M( N) =A LO G( RH M( KK ) ) +A LN R/ GM (K K) 

"rHmCn) =ExP (RhM (N > ) ■' 

G=2. «GM( KK )/ (GM( KK )-l, ) 

QM (N ) = SORT (QM( KK )»OM (KK) -G‘« ( PM ( N ) / RH M ( N ) -PM( KK )7RHM( KK ) )') “ 

WM (N ) = WM (KK ) + ( OM (KK) *QM( KK )-QM (N )«QM (N ) ) /2 . 

RM (N ) = RM (KK) ' " ' ' 

__ TM<M)=FT(PM(N) rRM(N) ,WM(M) ) 

GM (N ) = FGAM (TM( N) ,PM( l\l) ,RM( N) )■ 

EMM(N).= (QM(N)/SQRT (GM{N)<^PH(M )/RHM(N) ) ) 

XMUM (N ) = AS IN (1 »/EMM( N) ) 

Al = . 5^ (GM( N) /S IN (XMUM( N) )/C0S( XMUM (N ) ) +GM( KK )/SI N( XMUM (KK) )/COS 

KXMUM(KK))) ■ 

_ THM( N) =THM (KK) +0PT«ALNR/ A1 

1'2”cONTinUE '■ ■' ■ 

ER4 = THM( IFAN)-THM( IFaN + 1) 

IF(A0S(ERA),LT.1,E-O4) go to 15 ■ ' 

ITT=ITT+1 

IF ( 1 TT .GT. 15 ) GO TO 10 2 

IFdTT.GT,2) GQ_JO 14 

ER1=ER4 ■■■ 

BETlrBET 

BET=8ET+ .02 " 
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GO TO 3 


lO $ 
20 3 

wR iT (5( 6^ 2^3^ 
rORMAK* ERROR IN 

BETA LOOP IN C 01^ L« ) - 


STOP 


14 

DUM1=(3ET1-ER1«(BET 

-8ETl)/(eR4-ERl) 


ERHER4 ' ' 
BET1=BET 



BE T= DU Ml 
GO TO 3 

A 

15 

CONTINUE' 

RETURN 



END 






SUBROUTINE C0WL0<I'»THM, EM, BeTfPMMV 

COMMON/F VAR/ _____ „ 

i RH OF ( 2 00 > ,"CP XF ( 2 OO ) EM T ( 20 O') 7X MU F (' 20 0 ) , W F"( 2"0'0 ) , R F ( 2Q'0 > , G AM F ("20 0 ) 7 

2XMASSF (200|, HF( 2QQ) ,SF( 20 0) ,ALPP( 10 ) , 

3THETa(20) “ 

COMMON/XF/XF (200 ), YF (200 ), PF (2 00 ), OF (200 ), TF (2 00 ), THF( 20 0) ,ALPF( 

“ ■ HO ,2 00 ) 

CO MM ON /D / GA M I NF , E M I NF , R I N F , W_L!^£ 

qOMMoN/Ci'/'EmC I, TClV'OCl 

J=I+1 

'bEt=thf( i> »xmuf( D 

BET=BET- j. 1« ABS( 9ET) 

ITT=1 

3 CONI INUc 

cAlL Shock (Bet ,of(T) 7THr n )/gamf (T ), emfTT) ,Rhof( i) ,pf‘( i') ,wr( \u 

IRF (I ), QF (J ), TH F( J) ,G aHF(_J) ,EMF (J >, RHOF (J ), PF (J ),_WF (J )^F (J ), 

IJF (J), XMUF (J ), -1 .) 

P2 =P F< J) AM IN F« EM IN F« EM IN F _ ___ 

P2 =P 2/ PM M 
TH2= THF( J) 

'call PM(P2;iHM ;EM, TH2r,E'M3',~l) 

ER4==TH2T“TH2 _ _ 

IF (ABS (ER4 ),LT .1 .E-04) GO" TO 15 

1TT=1TT-»-1 _ _ _ __ 

IF(ITT.GT.15) go to’ 102 

IF(ITT,GT.2) go to 

'ER1 = ER4' 

BET1 = BET ^ __ 

BET=8ET+,02 

GO TO 3 _ _ _ 

102 WRITE(6,203) 

20 3 FORMAT (4 ERROR IN BETA LOOP IN COWLO<^) 

" ST OP • ' 

14 DU Ml =8 ET 1- ER (B ET -B ET 1) /( ER 4- ER 1) 

" ER 1= ER 4 ■ ■ ■ ^ ' 

BET1=BET _ 

'”bet=dumi 

GO TO 3 . 

"if cUnt’iNuE^^ 

DO 3 00 JJ=1,7 

300 ALPF (JJ, J)=ALPF( Jj,l ) 

xr ( j ) = xf ( i ) 

■ YF (J )= YF ( I ) ■ “ 

EMC1=EM3 

DU Ml = 1 , + ( G AM I N F- iT ) / 27 M I NF7 e'M I NF " 

DUM2=1 ,+(GAMINF-l. )/2, »EM3«EM3 __ _ _ __ _ ______ 

TC1=DUM1/DUM2 - - - - - 

QC1 = EM5/EM INF<^SQRT(TC1 ) 

RETURN ‘ ■" 

END 




— 2 » P2 P.ill hlA* U2>f^SH ) 

itVi 

_ _ VT=Q1«C0S( BET-THl) 

Ui=ABS (Ql^^SlN(BET-THi) ) 

__ GM1=G1-1 . : 

GP1=G1+1 . ' 

_XM 1 = El^il S I N ( BE J-lIll ) 

XMS=RH1«U1 

I F' ( I T , E Q . 1 ) U M = U 1« ^GM 1 « X M X M_1 + 2 . J / GP 1 .^X M 1 / XhU 

”■ 7 ~ RH 2= XM S/ UM"“ 

P2=XMS^^( Ul-'UM) +P 1 

VI =Ql<^Ql ' ' “ ‘ 

V2 =V To VT +UMOUM 

V V2T/'2. 

PH 12 =P Nil __ 

RH 2T =R HE Q ( H2 , P 2 r pH 1 2 7T 2‘)‘ ' “ ' 

ER=RH2-RH2T 

IP (ArS (ER/RHl) ;L'T. i. E- 03 ) "G0‘ TO” 9' 

IT=IT+1 

IF Tl T7GT .15)' G 0"T0 100 

IF(1T.GT,2) GO TO 11 ^ 

ER 2='ER^ 

U2=UM _ _ 

UM =0m«".9 9" ” " ' ' 

GO T0_7 

1 0 0“ WR' I T E ( 6T20 O' )~ 

200 format (* ERROR IN HUGON I OT_ LOOP IN SHO^KO 
■ STOP 

11 DUM2=U2-ER2« (UM-U2 )_/ (ER”ER2) 

ER2“ER \ 

U2=UM 

ijM =0 uM 2 
GO TO 7 

" " 9 "'IF (F SH ;G' T7"07 )' TH 2= 8E T- AT AN (U m7 VT ) 

IF (F SH ,L T. 0. ) TH 2= AT an <U H/ VT )+ BF^T 

Q2 =S QR T( V2 )" " 

G2=PGAM( T2 ,P2. PH 12 ) 

” ““ EH /S'OR 'TTgP oP 2/ RH 2T 

XHU2=ATAN( 1. /SQRT( EM2*EM2“10)_ 

^“"■rETURN^ ' 


SUBROUTI iME ShOcPK IS ,B£T',BeTN» A» B) 

COMMON/V aR/RHO{200 ), _ „ 

I'eM (2 o6)7xMu( 20 0) ;CPX (2O0‘rr‘W( 20 0) ,R(20O )VgAM( 20 0 )7X MA SS"(2 00 7 

2 XN (2 00 )» YN (2 00 >1 QN <2 00 ) , TN (2 00 )* PN (2 00 ), THN( 20 0) ,RHON 

3 (200 )7eMN( 20 0) rXMUN( 20 0) ,CPXn( 20o)',WN( 20 0)7RN( 20 0> ,GAMN( 20 0) , 

4 XMASSN (2 00 ), ALPN (1 0* 20 0) ,S 1( 10 )» HI (10) . _ ER R( 20 ) , TE HP (2 0 ) 

■■5 ,AlPdUM( lO ) ■ 

__ COMMON/X/ X( 20 0) ,Y (2 00 »Q^2 00 ) , T( 20 0 ) , T H( 20 0 )_r A LP (1 0, 20 

qO mh om /^- tx /X j; xJ 1 . mp tS , I o» iR ep l. ic he m; IP u7ch ;T de sg nV ir ,n XX ;x 8P‘. 

1 YBPf THBP ,R AD ,XBOD» YBOD »7 hBOO ,R ADB, XE ND ,T mend ,R TH ,YEX IT ,T HS T, TEST , 

1 iR pL ,Y 0» RA d8 2^ RR aD (2 0) f NS tA R* yN OZ ,E IN , PENrHl6,Hl7 

IP (6 ,EO. 0. ) BE TN=BET + ( THN( IS -1 )+XMUN (I S-1) -T H( H-1) -X HU (H-1 )) 

■■ ■ 

J= IS +1 

lms - ' -- " 

ITT=1 

3 LL S hO C>< (B El n 7 QN ( I ) f TH N( 1 > ,G aH (^( j ) , E MN ( I ) » RH ON ( I ) , PN ( I ) » WN { f ) , ” 
_ IRN (I ), QN (J )f THN( J) ,GAMN( J) ,E MN (J ), RH ON (J_), PN (J )r WN (J )^RN (J )_, TN (J ), 
IXHUN (J ), 1. ) 

EM1=XM1( .!5^.5.THN( J) ,XMUN( J> ,TH(M) ,XMU(H))__ 

■■”EH2=XM2(\5 7. 5. TH(M)EXHU(MyVTH(M+l> ,XHU(H+i))''' 

YDUM = (EM1<^Y( M) -EM2^»YN(L) ^EHlflEH2*( XN (L )-X(y^) ) ) / ( EMI- EM 2) 

RAT=(YDUM-Y(M) )/(Y(M + l)-Y(H) ) 

Y2 = YDUM _ _ _ _ _ _ 

IF(IS.E0,2,AND,M.E0.2) RAT = 0/ " ' ' 

IP(IS.EQ.2,AND.M.E0.2) Y2'Y(M) _ 

X ’ 2 = X rM)+RAT«()( (M+iy-X (MlT ' 

CALL PANDW (PJM ) , P (M + 1 ) ,_RAT , P2» DUMl, DUM2. DUM3 , DUM4, DUH5 , D_UM6» DUM7 

1,1) 

TH 2 = TH (M)+RAT«(TH (M + D-TH (M)) 

XMU2=XMU (M )oRAT# (XMU (M + 1 )-XMU( M) ) 

GAM2=G am (M )+RA T« (G AM (M-H )-GA M( H) ) _ 

' Al =P 1^ aT bV 'xM (j2 ,G aM 27 P2‘,X mU N( j) >G aM>k 7) ,P7‘irj) ) ' " ‘ 

A2 -f 2( A, R, 1. ,X J, XJ 1, X2 ,Y 2, TH 2, XM U2 ,X N( J) ,Y N( J) ,T HN (J ), XH UN (J )) 

' DUM1=A2*(XN(J)-X2) ■ ■ - - - 

DUM= (TH2-THN (J )-DUM1 )/Al _ 

PTEST=P2^^EXP (DUM ) ' “ ' ' " ' ^ 

ER4= (PTEST-PN( J) )/P( M) 

‘If (Ags (E'r 4 yVLT ;i ,E-0’4) Go''T0*’i5 “ 

ITT=ITT+1 

IF7I IT ,GTM5 )”G0 TO 102‘^”-' “ ‘ 

IFdTT .GT,2) GO TO 14 

"■ ER1=^ER4 ■ ' 

__ BET1=BETN 

' ■ BETN=BETN4 .02 ^ ■ - 

GO TO 3 . 

10 2 WRITE(6, 203) ’ ■ “ 

203 PORMAT(i^ ERROR IN BETA LOOP IN SHOCPT^^) 

STOP ” ■“ -- 

14 DUM1=BET1-ER1M8ETN-BET1)/(ER4-ER1 ) 

"■■ER1=ER4 ‘ 

BET1=8£TN 

b£tn=duhi ■ 

GO TO 3 

15 CO NT IN Lie " ■ ■■ 

__ XN(J) = XN(I) 

■ 'yN ( J ) = YN ( I )■ ' " 

DO 300 JJ=1,7 

30 0 ALPN (J J, J) =AlPN( JJ ,I)“ ' 
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‘ SUBROU TI NE P ANDW (P Ir P2 ,R aT ,P 0UT; wi".W 2# Q1 tQ2f Q3 rW QU T/IE Q. 1) 

DUM1=AL0G(P1) _ ; 

■" DUM2 =AL0G ('P2 ) 

_ POUT =DUMl+RAHMDU_M2;;^DUhl_) 

pO UT" =E XP (P OU T) 

IF (I ,EO, 1) RETURN 

WT 1= W-1 +0 Ql /2 . 

WT 2= W2 +Q 2^ Q2 /2 . : 

' wT 3= wT 1-^ RA T« ( W T2 -W T1 ) 

WOUT =Wi+RAT« (W2-W1 ) 

“ IF (I eQ .Eq. i> N0 Ut= wT3*'03ttQ3/2’- 
RE TU RN 

end 



Sii BR OU T I ME 'D R f E ST ( XN rV N , EMI , BE t » I D RO P i L" » M ) 

COMMON/X/ X(2 Q 0 ) > VJ 2 0 0_ ) » PJ..2 0 .0 ) > Q ( 2 0 0_)_» T ( gQ Q ) > T H ( 20 Q ) 

~EM"2= TA nTbE T> 

XN J= ( Y ( H ) - Y ( M+ L) -E Ml »X ( M )_♦ E M 2 « 

'yNT=’y<"M>'+£ Ml -ftTxM T“ Xt M> ") 

DI SNqR =SQR T( <XIM~X( M) )««2+( YfM-Y (M )) 2> _ _ 

D! SE X =SQRT( (X nT -X (M )) 2^^ (YnT ~Y (M )) *^^2) 

XF ( ( p I S^X^D r SN OR )_/_D I SE_jX* JJ 1 ) I DR O P =1 

“return ^ 


-A3?" 



SUdRoUtI K'P’ P- xH l» TH2» XM2» is ) " 

CO MM ON /D / GA M I NF , g Hj^NF 

G= 6 A M I 

GM1=G-1. 

GP 1=G+ !• 

DUM=1. +GM1/2 .-^XMlttXMl 

P=PiMt( GMi/G) ' *’ ' ' " 

_XM 2= ( D UM /P -1 , ) «2 Gf^l i 

XM2=S'0RT ( X M2 ) 

GX=SQRT( GPl/GMl) ^ - 

DUm1=SqRT( XM l^^XMl-i* > 

DUM2=SQRT( XM2«XM2“1. > 

DTH= GX<»( at an (DUM2/GX )- AT aN (DUMl/ GX ) ) +A TA N( DU Ml')- AT AN { DuM 2) 
TH2=THl+Fy)A_n IS )^DL)i 

■return 

END_ 





“ SUBROUTI NE' At_L ■( CPX. W^RHO, R, QAMVEM /XMU;q~t ;P /ALP V 
__ CO MM ON / 1 E_0 / 1 tu , P 1 N » R HO I N K , U1 NF #P I . 

~'common/a 7 TlN»'CPiN;RO 

CO MM ON /B/ GAMI NF ,EMI NF »R_IN_F, W1 NF 

common/B/"’wTmOlE - ■ - • • - 

COMMON /SP/ NS P 

COMMON VENTH/ HX 

DI MENS ION HI ( 1 0 )_,C PIJIJ ) .,D CP 1 (JLO ), ALPl ( 10) 

di'mEns ion' wt mo'Le (I'o”) ' 

CP X” 0 • 

iF'd EO .EQ'.'i'V G"b~"T'0 69 

_ W=0. 

" '"hx ^o , 

CA LL_T HE rM 0 ( T, ^:UC PI , ^ 

qO i 0 1 » nS p 

CP X= CP X+ AL PI (j )* CP 1( J) 

X-^ pi (0 )» Hi " ’ 

10 w=W+ ALPl (J )/wT MOLE (J)_ 

w= 1. /W 

RH0=GA^UNF*E_MINF»«2J|^ojVT/WI^NF 

■r=ro/w‘ 

GAM=CPX/ (CPX-R/CPIN) _ _ 

EM=Q«EHI NFttSQRt( GA M 1 NF «R‘IN F/ G A M/ R/Y) 

XMU=ATAN(1,/SQRT(E 2 " J 

RETURN 

69 J?HO = RHEO(W,PpR, T) 

gam-fgahTt,p,r) 

__ EM = Q/SQRT ( GAM«P/RM0) __ 

'XMU=ATAN(lV/SQRT(EMi»EM>-l', ) ) 

_ RETURN _ _ ^ „ 

END 


SUbRoUTI NE The RM0( T1 ,H ,Cp, DCp’) 

COMMON /SP/ MSP 

’ '"common/a7 tTn»cPin,Ro 
CO MM ON /B/ WTMOLE 

dimension WTmOLE(IO) 

dimension H(10 )f CP (10) ;DCP(10) __ _ 

T=TUT IN 

Cl=RO/CPIN ^ 

""c2=Ci/TIN 

_ _C3=C1^»T1N _ 1 

N= NS P 

DO 10 J=i»N _ 

H1=C2/WTm0L.E(J) " 

H2 =C 1/ WT MOLE (J ) 

*H3=C3/WTH0LE ( J ) 

CALL COEFF (J ,T ,A ,B ,C ,D ,E ,F ,G ) 

H( J) =T«( A+T* (B/2 .-►TttCC/G (D/4 ,^E/5, «T ) )') ) '*'F' 

H(J)=H(J)«H1 

■ CP(J)=A-«-T»(B+T«(C+T«(D+EeT)) ) 

CP (J )=CP (J )«H2 _ 

pC p( J) -B +T #( 2« ftC +T 3* *1^ +4 T) ) 

DCP( J) CP (J )« H3 
1 0 CP nT” 

RE TU RN 

^ ■ END ' . - - ^ 
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■ ■ ru NC T T ON ■ F T ( PT/F ^ H 5 ) 

COMMON/ I PP/IPP 

■ CO MM oN ‘ / TH E/ A 1 , a2 , A 3 » , A 5 » A6 I X MM 1 

CO MM ON / 1 EQ / 1 EQ » P I N » R HO INF> U i NF > P INF 
Cf^ MH oN /A7‘ Tl nVcP IN ,R O 

pATA 163/0/ 

iF lA G= 0 

iHpL_D= 0 

p= Pi /p iN iN po i , bl 32 +0 5 

H=H5«U INF/10^7 63 9/ 1. E/06 __ 

'pS^F^F " " 

A10=ALOG(P)/2.3-5. 

"29=. i2 5«a10<^a 10' ‘ ".2 75<vAiO 

IT =1 

l"F (T 63 ,E"q. IT'G 0 fO i 60 0 

__I63=1 _ 

T=15oO. 

TO =1 50 0^ 

'“"’if (F .GE. 0. ) GO TO 12 0 

T=60 0. 

fO = T 


100 0 CONTINUE _ 

IF (F '.LT. 0 . ) GO TO 400 

GO TO 120 

' 50 * E 0 = ( H-H 1 )/H‘ " ' 

__ I F ( A B S_( E 0 _) . L T . I . E- 04 ) G O TO 340 

"Too ’t"~=toci 7 i 

__502 IT = 2 _ 

IF (F .LT . 0 . ) GO TO 400 

__ _ GO TO 12 0 

lO"o El=(H-Hl)/H 

IF_(A BS (E IJ ,L T_^ 1 ._E- 14 J 

“"lT=ft+f~ 

IF (I T. LT . 21 ) GO TO 10 _ . 

IF (AbS (T -2 00 0 . ). LT .1 0 . ) GO Tq 63 0 

IF (I HOLD ,EQ. 0 ) GO TO 8 00 _ _ 

■’ WRITE! 6 i 831 ) Pl»H 5 VT” 

831 format ERROR I N F T «/ P I = ^ HI 

i'T 

STOP 

80 0 ” T=TH 0 LD 

H=,HH 0 Lb”~ 


«El3.5»5Xf^T = -»E13.5 


IT = 1 

IHOLD=l 

IF(T.L*E.2000. ) IHOLD = -l_ 

GO' TO 1000 

83 0 IF ( I FL AG . E 1/ G_^ T .? 

iflag=i 

_ T0 =2 00 0_^ 

. - 20 oO . 

IF (F .LT . 0, ) ^GO TO 40 0 

GO T 0 12 0 " ^ 

50 4 WR IT E( 6j 11 ) El 
11 FORMAT (» temperature IN FT SET 

GO TO 340_ 

IQ T9=T0-E0*(T -Td)/(E1"E0) 

505 E0=E1 . 

to=t 


TO 2000’ ’ ^ERROR '=’ <T‘13’, 5) 


AA3“ 



T= T9 

If (f ,LT, 0. > gO Tc 400 

1^2 0 A= E- 07 «( 10 42 2 24 2* F+ . 9 67 ) 

’ B= . 0 01 . 0 11 67 »F 2 +♦ 15 03#F + , 93 8) 

C=-» 0264i^F2 +. 6731 »F + . 4293 

■ IF (F .LE. 1. ) GO TO 190 ‘ 

A=l. E-07»< 1. 78 7«F2 - 5. 48 ^»F -»-5 . 4 ) 

B= .001*( -.1867«F2 +l.ll»F+a76) 

C=-.0933<jF 2 +3 .975*F-2 ,8Q8 
'l9‘0' IF (t ,LE, 20 00 . . AND. IHOLD. NE ,1 V GO " TO ‘29~0' 

A= .0Q0 001»(1 .792«F2 + , 39 03 «F + . 31 ) 

8=^ .0 01 “9 ,0 F2 - .0 7917 +. 24 5) 

C=10 ,86«F2 - .1163*F<^ .97 
■ IF (F .LE. 1. ) GO TO 290 

__ A= . 000001» (4 .81«F2 - 13 .9 «F +1 1 . 59 ) 

B-’iO 01 -23. 08 «F 2 -^6 6 . 82 »F'-5 2/61 ) ' 

C=27,05»F2 -73 ,73i»F+58 .39 
29 0 HI =A -»T cT +B «T -j-C 

IF (T .LE. 20 00 .. AND. IHOLD. NE ,1 ) GO TO 370 
HI =Hl» (1 . + (1 ,^-F) tf ( T/20 00 ,-l, )«Z9 ) 

37 0 CO NT INUE 
GO To 35 0 

40 0 T2=Ti^T _ _ _ 

T3 rT 2« T " " " ‘ 

T4=T3«T 

T5 =T 4^ T ' ' 

IF (F .LT. -1.5) GO TO 45 0 
■ XHMT=16'.'043 “ ' " ■ ■ 

A1 =4. 2497670 
A2=-6. 9126562E-03 
A3 = 3 . 1602134E-05 

A4=-2. 97 15 43 2E -0 8 

A5 =9 .5 10 35 B0E-12 

"'• a6 =- 1- o1’86'63 2? +0 4 

GO To 46 Q 

45 0 CO NT INuE ‘ 

A1 =1 .1 20 24 36 

A2=1 . 3 90 57 16 E- 02 

A3 =2 .6 56 83 74 R- 06 
"A4=-i. 15 60 27 2E -08 ^ “ 

A5=5 .2386929E-12 

A6 =5 .3328896E+03 

XMM1 =28. 054 

46 0 H1=A1« T+ A2»T 2/2, +A3«T3/3 . + A4 4/ 4 . +A 5» T5 +A 6 

H1=H1^j 8314 ,/XMMl/l .E+06 _ _ 

35 0 IF ( I T. EO .1 )■ GO TO 50' - - 

GO TO 100 

34 0 T0=T ■ ■ ~ ■ 

FT=T/TIN 

THOLd=T ■ ■ 

hhold=h 

■"'return . 

end 
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rUNCTION'T'Hi 

COMMON /THE/ A 1 » A2 » A3 r_A4 , Af5 , A6 _i XMMl 
CO MM ON /rEQ'/rEQ' r P'Wl HO i N rruT'NF ; P'i N F 

CO MM ON /A/ TINiCliINjRO 

p-pl /P jN *P IN i. bi 32 5E +0 5 


T=T 1 IN _ 

F 2 =F 

(P_.L T. 0 ^ )_^ _GO_T 0_40 0 __ " 1 

iFTr;GTT 2 oob , ) go to 190 

IF (F ,GT, 1 . ) GO TO 19 1 _ 
i 2 l '' a= E- 0 '/ -'♦’iO 42 <»F 2" + .« 24 2 <^F+ ,987 ) 

_ 8 = ,0 01 ,0 1167 «F 2 +• 1 ^? 0 ^^F + , 93 8 ) 

Q= -. 02 84 «F 2 +, 67 3 I ^F +. 42 93 

G 0 _ TJ ^_29 0 __ _ 

l 9 1 'a= f. E’ 07 <rri\ 78 7<^ F 2 - 5 . 4 B *5 J ) 

B= .0 01 i^( 18 67 *F 2 +1 . 1 1 ^ F+ .1 76 ) 

■ c= -- 09 33 *F 2 +3 .9 75 ^F -2 ,8 06 

GO TO 29 0 

19 0 ‘ IF'(F ;GT. 1 . ■) GO TO 19 2 

A= .0 00 00 1 ^ (1 .7 92 <*F 2 39 83 «_F_+. 31 )_ 

^ 9=70 Oi » { - 9 ; 0 F 2 “ - 70 79 i' 7 “ttF + ; 24 " 5 ) "■ 

C= 10 . 86 «F 2 -_. 1183 ^^F+ . 9_7 

“ GO 'TO^ 29 0 

192 A= . U 00001 « ( 4 . 81 «F 2 -13 » 9 »F + 11 . 59 
■ B = . 001 ft( “ 23 , 0 O»F 2 + 66 . 82 - 5 ^F -52 .' 6 i >' 

C = 2 7 , ^5 ^ F 2 j- 7 3 . 73 ff F + 56.39 

290 H 1 = A » T T + B t + C 


IF( T . LE .2000. ) GO TO 370 

■‘A10 = ALOG(P)/2.3^5 7 
79= . 125« A10*A10 _ -,_275«A10 

H1=H1« (1 . + (1 , + F) «rT/2b06 ,-i7 )»29 ) 

37_0„ HI Ij 

GO ^^ b 34 0 


40 0 T2 =T ftT 
T3'=T 2i» T 
T4 =T 3^ T 


t5 =T T 

HI =A T+ A2 <>T 2/ 2 ._+A3<^T3/3 . + A4 ^ T4/ 4 *, +A^5^^A 6_ 

H1"=H 1» b3'i 4 , / xM Ml 

34 0 CONT INUE 

' 'FH^Hl-alO .7 63 9/Ul NF 

RETURN_^ 

END 

function FGAM( T1 ,_P1, F)^ _ 

■■■■“■ ■ COMMON /THI/' aI, A2 ,A3» A4 A6 ,XMM1 

COMMON/IEQ/IEQ ,P IN rRHO INF, UlNF ,_P INF 

"comm ON /A/ T I Nr CP IN rRO 


__ _t=ti«t in 

T2=T «T ■ ■ 

P = PI /P I N I NF « 1 . 01 32 5 E +0 3 

XMTTbT ' 

IF (F .LI. 0 . ) GO TO 550 

IF (T ,LE. 10 00 . )■ GO TO 4 4'Ci" 

XM=-2. 15E“0e<»T 2 +. 0 0 0091«T-, 0695 
44 0 XN’=4 .E-0 9«T2 - . 0 00 02 ifT - . Cl 9 

IF (F .LE, 1. ) GO TO 47 0 

XN"=7q3 39 «-S qR T ( T )'-7 0O'o 3 91 i^T 68 1 

470 G=-l .833E-07«T2 +. 00 00 75 *T +1 ,3 67 
IF (T .LT. 500, ) GO TO 52 0 

-aA 5 " 



G=2.E-08«T2 00 0133»T^-1 .423 

“ ■ IF (T X^• 20 00 .) Go TO 520 
_ G=7,267E'-08<»T2 - . 0 00 45 7«_T<- 1 , 8 

52 d” G= XM aCoG (P')/ 2. 3"5*’ )* XN iCf 

GO TO 53 0 
55'0 T3=T2»T 
T4=T3«T 

"CP=A1 + A2-5^T+A3<^T2 + A4<j-T3 + A5«-T4 

_g=cp«<j:P“1. ) ^ 

53 o’ CONTINUE ■ - - 

FGAM = G _ _ ' 

" "return “ 

END 





rUNC l'l ON' RHIt'qY H/P1 ff' »T ) 

COMM 0 N / i BQ / 1 tU , H I N # H HO I N r , 'J I NF >P I NF 

COmKoK' /AV" TI'nVcP IN ,R 0 

T1=FT(P1^F_,_H_) 

' T="fr«t”iN 

P=:P1 /P 1N«P INF^^ 1. (U32 5E+0 5 

iF (F ,L T» 0- ) GO ”^0 22 60 

FN M= 1 . 53 -b-F *F -5 .8 95 *F +2 8. 

'pN N= Y»'6* F~*T~ lO .6 «F +3 3* 6 

I F ( T , G T . 20 00 . ) G 0_ TO __2 03 0_ 

XM =F hM 

JF(F,LT.1._) GO_T_q^_2160 ^ _ 

■■■"■xm=Fnn “ 

__ GO^ T 0 2 1 6 0 

2 03 o" fF =F *'F 

A=~2 .3«FF+4, 01«F+1 ,7 36 ____ 

B= 8-"6l F- l5 .4 2*^ F' 6* 66 

C=-16. B8ttFF+33 ,2 1^»F+14 ,5 8 _ 

”XN =- .4 37 5*^ FF + . 06 25Y»F +2 .0 8 ’ 

D = A^^ ( ALOG (P) /2 .3 .5 + B^t- ( ALQG ( P ) / 2 . 3 ) ->-C 

XM =F NFr^D Vn T -2 00 0’» )/ iO dO“'.T«« XN 

IF(F.LT,1.) GO TC_ 2160 

A=-. 62 2*FF+2 .3 63^F+1'.9 05' 

B = 2.76«FF-*7,56«F-8.66 

C = “3'.6wFF + 7.36<^F+27.l5' 

XN = - .4;MFF + 1 .825tfF + ,_35 

D = A»( AL0'G<P )/2 . 3) »»i. 5 + B« ( ALOGUM /2 . 3 )“^C 

X M = F N N - D ( T - 2 0 0 0 . ) / 1 0 0 0 , ) » « X N 

"go TO 2160" 

_2260 kF=F-.5 , _ 

IF(KF.eO.-l)XM=16. 043 

IF (KF. EQ .-2) XM =2 8. 05 4 

2 i6 o' rH pQ rP *X fV t/83'i4'.3 »6 '.2'42 GE -0 2/ 32 .1 74 

T= T/ TI N 

RHEQ'=RHEQ/RHOi NF 

^RETURN ^ 

end" 



rUNC Tl ON XMl (ALPHA ,8 eta* TA ,Xa, TC ,XC) ' 
XM1= alphas TAN ( TA +X A) 

If (BET A« GT .o', ) XMl=XMi+8ETA*T aN (T'c+XC ) 
RETURN 

end"" 



rUfC- Tl ON XM2 (Ai;n, TA ,X A; TC ,Xc) 

XM2= AL «T an (T A-XA ) __ 

*iF'(B7G'T . 0 » ) X M2 =X' M2 +8 AN ( T C- XC ) 
RETURN^ 

end' 



FUNCTl ON Xm 3 (A ,B ,TD^ fC ) 

XM3= A«TAN( TD )_ 

IF'(B ,G T /O .'0 > XM 3= XM'3+ 8» tAN( TC ) ' 

RETURN 

END 




fVnC iT 0N"F i'(A* B'," XHuiV GAMi , 'pi ;xm'U2 , 

lGAH2i P2 ) 

Vi =STn '( X MUl ) *C OS ( X MU 1 j /G AM i 

IF (B .G T. 0. 0) FI =( FI +S IN (XMU 2) 

ICOS< XMU2 )/GAH2 '>72. 

RETURN _ - 

■ END 



FUNC TI on ^ 2( A» 8» OP T# XJ ,X Jl ,X ,Y ,T"h, XMU» XN ,Yn, TH N, XM UN ) 

F21=0. _ _ 

ipIxj.fo.o.) Go to 15 

F21= A<»SI N( TH )»SI N( XMU) /cos {Th+OPT»XMU) /Y 

IF (B .GT. 0- ) F2l=r2l‘"B*»Si N('tHN) «S IN (X MU N)'/C OS {ThN +O PT »XMUN) /Y N 
15 F22=0. 

iF (X jl ,£q. 0 . } GO TO lO ■■■■ ■ 

F2 2= A#COS( TH )«SI N( XMU) /COS (T H+ OP T« XMU) 

IF (B .Gfi 0.’ ) F2 2=F2 2+B^C0s( THn) »S JN (XMUN)'7’C0S (Th'N VO PT Vx MU nF 

10 F2=F21+F22 

rB TOrN ■ '■ ■ 

END _ _ _ _ _ 




_.SU BR ou T I NE X MASS S ( R A Tli .^1^1 S > 

COMMON /EVA R/' 

IRHOr (2C0 ) , CjPXF ('^OU )jJlEL2p_0_) ,XMUF( 2Q0) /Wf( 200J_.RE UriAiltilQll/. 

‘2XMAS"sFr2 0Tn » “ H F { 2Q 0 > , S F ( 20 0 ) , ^LP D ( i 0 ) » 

3THETA(20) . . .. 

COMMON/XF/XF (2 00 ) , YF (2 OC ) » PF (2 00 ) » QF ( 2 00 ) . TF (2 00 ) , TH F ( 20 0 ) »ALPF( 


110 f 2 00 ) _ _ 

do" 170 "l =l',NPTS 
C A_LL__F M_(_F HI ^E MF (I ) , GAHF (\ ) ) 
"fM2=FM 1/RA'TH' 


nM=i __ 
EMT=RA TH^^EMF ( 1 ) 


1 9 7 CONTINUE _ _ _ 

C X L L ' FlX( F M T , EM t ; G A M F” { I ) ) 

ER M=_( F H2j’F HT ) / 

TF ( A B S' ( 't R H") . L T . iT^-Q STTFO T 0 T7 1 

PUMD=1_./RATM _ _ 

" ffM= jTM+T ■ ' 

_[F U TM .GT. 15 ) GO T 0_2 

\f U TMVG‘t:2) go' to 1 

EM Tl =E MT 

ER Hi~^c RM 

EH T= EM TiM'JUMD ^ 

GO' To" 197" 

1 DU M= EM Tl -e Rii 1«_(E MT ~E MT 1) / (_ER M- EH Ml ) 

EHTI-EmT"' 

ER Ml =E RM 

emt^c-Um" 

GO TO 19 7 

2 WRITE(6>'5) 

3 format ERROR IN ITERATION LOOP IN XMASSS*) 


STOP 

1.71 continue 

’• "r.lJ QATTF (I )+l» )/ (G^riK^ [ ) -1 .) ) 

ET SQ =S OR T ( EM T<^ EM T- 1 . ) 

FS'Q= SO rt (E jr ■( ry^ifcM Fn )--i7> 

XN U2 =0 UM TA R( ET SQ /D UH )- AT AN (E TS Q) 

■ xN Ul =D IJMYt-A TA H( eS Q/ PU M) -A TA N< ES Q) 

0 M U = X N U 2_- X_N U 1 ^ 

^CALC 'P HI (On’G XI ) 

17 0 CONTINUE 

return'""" 

END , 




SUBROUTING rM(rMC*XH,G) ■ 

rMC=XH/( i . + ( G-1. )/2. ) » » { (G +1 , ) /2 , / ( G “I , ) > 

■return 

END 





SUBRoOtI NE 'P^ * 

Cp .MM ON /f VA 9/ 

^i'rhoPT 2 0 oT; cP X F"r2'uolTeN r(‘20'oT7x )Tw'r rfoTTTRTT^ olj'rTGTt^^ 

2XM ASSF (200 ), __ _ _ _ 20 0 ) , S F ( 20 0 ) , A].P D < 10 ) , 

3THETa('20 7" ■ ■ 

COMMOlV/XF/^(2 00_7 ^_(2C0 ) , PF (2 00 ) , QF (2 00 )j_lF (2 00 7 JHFp 20 0 ) ^ALfF^ 

1 10 ,2'00') ■ ' ■ ■ 

DTH=2./5 7.3 

IFaN=DNU/DTH ■ - - .. - - -- - 

IF ( D TH . G E . ON U ) DT H= ON U 

IF rp'Thr;EQ. DNpyiF ^ - - - - - - 

P1=PF(J6^_ _ 

” 'T1-TF‘( J6 ) ■ '■ ■■ ■■ 

EH 1= EM F ( J6 ) 

'hi-WT’(J61 ■ ■ - — - 

G1 =G AM F ( J6 ) 

RH 1= RH OR’3 6 ) 

PH 11 =RF(J6) 

"' Ul=0Frj6 ) ■■ ~ 

Ul=Ul«Ui 

i7r =Hi+ui /2 r ■' ■ * ■ ■■■ " ■■ ■■ ” ‘ ■ 

PI =ALOG( PI ) 

■"‘XNU=0. “ 

XNUl =0 . 

DO llT r= itif an 

IF (I ,E0. IFAN )DTH-D!NU~XNU 

RKui’^TsTFriTyEMi y~“~ 

__ B1 =G1-^EM 1/COS ( XHUl ) 

>2=-8i i^DTFT+Pi ■ - ■ 

RH2= (P2-P1 )/Gl 

rH2=RH1-EXP( RH2)“““ ' 

P2P=EXP(P2) 

■pi p= l^_X p( pi } ■ 

U2 ~U I- 2. 1/ (G 1“ 1. )* (P 2P /R H2 -P IP /R Hi ) 

H2T^h T- U272 7”~ 

T2=F T( P2P» PH II ,H2) 

G2 -F G a M ( T2 r P 2P # P H I 1 ) 

A2 -G 2-:» P2P/RI] 2 

7^’ 2~= SQ Rr^FTp? A 27" ' 

XNU= XNUl +DTH 

XN’U1=XNU ' 

P1=P2 

"Ti =^2’ ~ " “ ‘ 

G1=G2 

"EH‘i=EM2" ■■ """ . . . 

U1=U2 

RH1= RH2 ' ■' ■ " 

10 CONTINUE 

p2=bxP(P2>” " " ' 

U2 jS OR T ( U2 )_ 

__ T F (J6) =T 2 • 

£M V (_J 6) =E M2 

RH OF ( J 6 ) =R H2 

PE TU RN 

END ■ : 

-A55- 



